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medium has been provided with the ability to obtain a com- 
plex permittivity E or a complex permeability p of a sample in 
a cavity. One or more complex -valued resonance frequencies 
(f,„) of the cavity, wherein each f m is a measurement, are 
obtained. Maxwell’s equations are solved exactly for E, and I 
or p, using the f„, as known quantities, thereby obtaining the E 
and/or p of the sample. 
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TECHNIQUE FOR PERFORMING 
DIELECTRIC PROPERTY MEASUREMENTS 

AT MICROWAVE FREQUENCIES 

CROSS REFERENCE TO RELATED 5 

APPLICATIONS 

This application claims the benefit under 35 U.S.C. Section 
1 1 9(e) of the following co-pending and conmionly -assigned 
U.S. provisional patent application, which is incorporated by 10 
reference herein: 

Provisional Application Ser. No. 61/240,564, filed on Sep. 

8, 2009, by Martin B. Barmatz and Henry W. Jackson, entitled 
“NEW TECHNIQUE FOR PERFORMING DIELECTRIC 
PROPERTY MEASUREMENTS AT MICROWAVE FRE- 15 
QUENCIES.” 

STATEMENT REGARDING FEDERALLY 

SPONSORED RESEARCH AND DEVELOPMENT 

20 

The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-517 (35 U.S.C. Section 202) in 
which the Contractor has elected to retain title. 

25 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

This invention relates to a method for obtaining dielectric 
and magnetic constants for a sample, and an apparatus for 30 
performing the same. 

2. Description of the Related Art 

(Note: This application references a number of different 
publications as indicated throughout the specification by one 
or more author names and/or years within parentheses, e.g., 35 
(Name, year). A list of these different publications ordered 
according to these reference numbers can be found below in 
the section entitled “References.” Each of these publications 
is incorporated by reference herein.) 

Most microwave dielectric property measurements rely on 40 
a cavity perturbation approach that requires very small solid 
samples (Meng et al., 1 995). Cavity perturbation methods are 
only approximate and have several assumptions that have to 
be satisfied to obtain accurate results. In the conventional 
cavity perturbation approach, it is assumed that the sample is 45 
sufficiently small, so that the perturbation in the electric and 
magnetic fields with the sample inserted in the cavity is neg- 
ligible. It is also assumed that the difference between the 
quality factors of the cavity, due to the cavity wall loss with 
and without the sample, is negligible. These assumptions 50 
limit the sample volume to be only a few percent of the cavity 
volume in order to obtain accurate results. 

SUMMARY OF THE INVENTION 

55 

The present invention has performed the additional calcu- 
lations that enables one to arrive at an even more accurate 
determination of the complex dielectric constant (more accu- 
rate than a cavity perturbation approach, for example). This 
improved approach is particularly important when working 60 
with liquid samples, larger solid or liquid samples, or samples 
with laiger dielectric constants that further perturb the cavity 
electromagnetic fields. 

The present invention may solve Maxwell’s equations for 
dielectric and magnetic properties by knowing the cavity 65 
quality factor before and after the sample is inserted, and then 
matching the solutions of Maxwell’s equations with experi- 


2 

ment. The Maxwell’s equations, applied to this invention, 
may be solved using a unique root finder approach. 

Thus, to overcome the limitations in the prior art, and to 
overcome other limitations that will become apparent upon 
reading and understanding the present specification, the 
present invention discloses a method of obtaining at least one 
property selected from a complex-valued dielectric permit- 
tivity (E), and a complex valued magnetic permeability (p) of 
at least one sample, comprising obtaining one or more com- 
plex-valued resonance frequencies (f„,) of the cavity wherein 
each f m is a measurement, the real part (Re) of f„, being the 
measured cavity resonant frequency, and the imaginary part 
(Im) of f„, being related to the measured cavity quality factor 
(Q,„) by Im (f m )=Re (f„,)/2 Q m , and solving Maxwell’s equa- 
tions (Sphicopoulos et al., 1984), after including wall losses, 
to then solve exactly for E or p, using the f m as known 
quantities, thereby obtaining the E or p of the sample. 

Both £ and p are complex valued such that E=Re(E)+jIm 
(E) and p=Re(p)+jIm(p). 

The method may comprise obtaining one measured com- 
plex-valued resonance frequency f m =i cavity of the cavity when 
empty; obtaining one measured complex -valued resonance 
frequency f „=f sample of the cavity containing the sample, 
wherein i cavir> . and f samp!e are measurements obtained by 
exciting the cavity with one or more cavity resonant modes; 
and solving Maxwell’s equations exactly for the E or p, using 
the i cavity and the i sample . 

The exciting of cavity modes may comprise exciting cavity 
resonant frequencies to produce electric and magnetic fields 
in the cavity. 

The method may further comprise obtaining one measured 
complex-valued resonance frequency i holder of the cavity 
comprising the sample holder without the sample. 

The cavity and the sample may be annular or cylindrical 
regions. 

The solving step may further comprise solving the Max- 
well’s equations for a complex-valued resonance frequency f 
of the cavity including the sample, wherein the f may be a 
function of the E, and/or p. 

The solving step may further comprise solving the Max- 
well’s equations for a complex -valued resonance frequency g 
of the cavity including the sample holder without the sample, 
wherein the g is a function of the ^holder- 

The cavity may comprise one or more additional materials 
having known material dielectric permittivities and known 
material magnetic permeabilities, such that the sample’ s E or 
p are the only unknowns during the solving step. 

More generally, the method may further comprise solving 
Maxwell’s equations for complex-valued resonance fre- 
quency f({E„ p,}), where a complex valued dielectric permit- 
tivity E, and complex valued magnetic permeability p ; are 
assigned for each zone i in a set {i=l, 2, ... N} of annular 
regions and a central cylinder, the set of zones {i} comprising 
the interior of the cylindrical cavity (Sphicopoulos et. al., 
1984). 

The method may further comprise obtaining an experimen- 
tally determined electrical conductivity ct w of the cavity’s 
walls using the f cavity . The electrical conductivity a w of the 
cavity walls may be obtained using the f cavity and known 
formulas for the complex-valued f of the empty cavity (Wil- 
son et. al., 1949). 

The present invention may use the a w to modify the f and/or 
g to account for a power absorbed by the cavity’s walls, 
thereby obtaining a resultant complex valued resonant fre- 
quency F (from f) that is a function of the E and/or p, and/or 
obtaining a resultant complex -valued resonant frequency G 
(from g) that is a function of the E /;oW( ,,.. 
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More generally, in order to achieve the modification(s) for 
f({G ; , p,.}), the method may further comprise using the sur- 
face resistance formula evaluated using the calculated mag- 
netic field corresponding to the frequency f({G„ p,}), and the 
electrical conductivity ct w of the cavity walls to obtain the 5 
power absorbed and corresponding quality factor Q w; 
for the cavity walls, and then obtaining the contribution 
fw 7 =Re (f {G„ p,})/2 Q wo „ to the imaginary part of the 
complex-valued resonance frequency, so that the resultant 
complex-valued resonance frequency being F({G,-, p,})=f 10 

(fe p,})+j f 

The present invention may fiirther comprise obtaining the 
G or p by solving F=f sa /e . Alternatively, tlie method may 
further comprise obtaining the Ei holder by solving G=f hoIder 
and then solving of F=f JO le using tlie E. holder obtained by the 15 
solving of G=f holder . 

The method may further comprise obtaining a range of 
values for the G or p, and selecting the value in the range that 
satisfies F=f sa !e to any degree of precision desired. 

A sample holder may be needed when performing dielec- 20 
trie and/or magnetic measurements on liquids or granular 
materials. Generally, the sample holder is made of a nonmag- 
netic dielectric material, such as Teflon or quartz. In many 
cases, the dielectric properties of these materials can be 
obtained from the manufacturers. However, when that infor- 25 
mation is not available, the dielectric properties of the sample 
holder can be obtained using the present invention. More 
generally, the method may further comprise obtaining the 
complex-valued dielectric permittivity of an empty nonmag- 
netic sample holder G holder in a zone i=I, or in two non- 30 
adjacent zones I and K occupied by the holder, where the 
dielectric permittivity G, for each of the other zones in the set 
{i} is known, by solving F({G i })=f m , where f m is the com- 
plex-valued resonance frequency for the cavity containing an 
empty sample holder. 35 

The method may further include obtaining the complex- 
valued dielectric permittivity E. sample of a nomnagnetic 
sample, in a zone i=J, contained in a sample holder, where the 
complex- valued permittivity G, for each of the other zones in 
the set {i} is known, by solving F({G ( })=f m , where i m is the 40 
complex-valued resonance frequency for the cavity contain- 
ing a sample in a sample holder. 

The sample may be large and/or lossy enough to perturb 
cavity modes of the cavity. For example, the sample may fill 
between 1 % and 95% of the volume of the cavity. However, in 45 
the present invention, the sample may be of any size relative 
to the cavity volume. 

Transverse Magnetic (TM 0 „ 0 ) modes, wherein n is a posi- 
tive integer, may be used to excite the cylindrical cavity. The 
integer n may be, but is not limited to, between one and three. 50 

The sample may be non-magnetic, have any size, and be 
positioned along the cavity’s axis. 

To obtain the p of the sample that is magnetic, the cavity 
resonant modes may comprise at least one TM 0 „ 0 mode, and 
a magnetic field and an electric field associated with the 55 
TM 0m0 mode. The method may then fiirther comprise per- 
forming the steps of the present invention for the sample that 
is a first cylindrical sample positioned along the cavity’s axis, 
the first cylindrical sample having a first radius that is suffi- 
ciently small so that the magnetic field’s strength at the sur- 60 
face at the first radius’ end, is negligible, thereby allowing the 
Maxwell’s equations and the F to be solved for only one 
unknown, the dielectric permittivity of the sample ( G sampIe ). 

Next, the steps of the present invention may be performed 
for the same sample that is a second cylinder comprising a 65 
same material as the first cylinder, and having a larger radius 
than the first radius, such that the second cylinder positioned 
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along the axis experiences both the magnetic field and the 
electric field, thereby allowing the Maxwell's equations and 
the F to be solved for only one unknown, the p, using the G 
obtained with the first cylinder. 

If it is desired to use only one sample for measuring both 
the permittivity and permeability, the present invention may 
place an annular sample (or a liquid or powder in an annular 
holder) located at a radius value associated with the maxi- 
mum electric field for the TMg 20 mode (where the magnetic 
field is zero). Measurements for this TM 020 mode may deter- 
mine the sample permittivity G ja 7e . Then, using the TM 0 , 0 
mode, where the same sample may be in both the magnetic 
and electric fields, the present invention may back out the 
sample permeability knowing its permittivity. With this 
approach, the present invention may assume that the permit- 
tivity does not change significantly between the frequencies 
of the TM 020 and TM 010 modes. 

The method may further comprise using the G or p to 
interpret remote sensing data obtained from extra-terrestrial 
or terrestrial material, thereby identifying one or more com- 
positions, optical properties, or dimensions of the material; 
and transmitting, storing, or displaying the compositions, 
optical properties, or dimensions of the material. 

The solving of F=f sampIe further may further comprise find- 
ing roots of F-f, 7e =0, by: (i) specifying a range of values 
for z within a first rectangle in a complex plane, wherein the 
z are complex numbers in a complex plane; (ii) introducing a 
mesh into the first rectangle, thereby forming the mesh com- 
prising mesh points at comers of second rectangles, wherein 
the second rectangles are smaller than the first rectangles; (iii) 
evaluating Re F(z) at each of the mesh points; (iv) determin- 
ing one or more first sides of the second rectangles where the 
Re F(z) changes sign; (v) evaluating Im F(z) at each of the 
mesh points; (vi) determining one or more second sides of the 
second rectangles where the Im F(G) changes sign; (vii) 
selecting the second rectangles having both the first sides and 
second sides, where both the Re F(z) and the ImF(z) change 
sign, to obtain selected second rectangles; (viii) interpolating 
the Re F(z) in the selected second rectangles to define a Re 
F(z) null line; (ix) interpolating the Im F(z) in the selected 
second rectangles to define an Im F(z) null line; (x) finding a 
first intersection point between the Re F(z) null line and the 
Im F(z) null line; (xi) setting the second rectangle as the first 
rectangle and repeating steps (i)-(x) using the second rect- 
angle as the first rectangle, thereby finding a second intersec- 
tion point; (xii) obtaining a difference between the first inter- 
section point and the second intersection point, and 
comparing the difference with a precision parameter that is 
equal to a desired precision; (xiii) repeating the method until 
the difference is at least as small as the precision parameter, 
thereby obtaining a final intersection point; (xiv) and setting 
the final intersection point equal to the property. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Referring now to the drawings in which like reference 
numbers represent corresponding parts throughout: 

FIG. 1(a) illustrates a simple embodiment comprising a 
cylindrical sample along the cylindrical cavity’s axis, accord- 
ing to one or more embodiments of the present invention; 

FIG. 1(b) is a cross-sectional schematic of the cylindrical 
sample in the cylindrical cavity of FIG. 1(a), according to one 
or more embodiments of the present invention; 

FIG. 1(c) is a cross-sectional schematic of a larger cylin- 
drical sample in the cylindrical cavity of FIG. 1(a), according 
to one or more embodiments of the present invention; 
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FIG. 1(d) is a cross-sectional schematic of the cylindrical 
cavity comprising an empty cylindrical sample holder, 
according to one or more embodiments of the present inven- 
tion; 

FIG. 2(a) illustrates a cylindrical cavity containing an 5 
annular sample, according to one or more embodiments of the 
present invention; 

FIG. 2(b) is a cross-sectional schematic of the cylindrical 
cavity of FIG. 2(a) containing an empty annulus holder, 
according to one or more embodiments of the present inven- 10 
tion; 

FIG. 2(c) is a cross-sectional schematic of the cylindrical 
cavity of FIG. 2(a) containing the annular sample in the 
cylindrical annulus holder; according to one or more embodi- I5 
ments of the present invention; 

FIG. 3 is a cross-sectional schematic of a cylindrical cavity 
for the theoretical approach for calculating the complex reso- 
nant frequency (Sphicopoulos et. al., 1 984), that may be used 
in one or more embodiments of the present invention; 20 

FIG. 4 is a flowchart illustrating a method of obtaining a 
complex permittivity E or complex permeability u of a 
sample in a cavity, according to one or more embodiments of 
the present invention; 

FIG. 5 is a flowchart illustrating a method of obtaining a 25 
complex permittivity E or complex permeability p of a mag- 
netic sample, according to one or more embodiments of the 
present invention; 

FIG. 6 is a flowchart illustrating a rootfinder method 
according to one or more embodiments of the present inven- 30 
tion. 

FIG. 7 is a schematic illustrating a mesh on the complex 
plane, according to one or more embodiments of the root- 
finder method illustrated in FIG. 6; 

FIG. 8 is a cross-sectional drawing of a microwave cavity, 35 
according to one or more embodiments of the present inven- 
tion; 

FIG. 9 is an exemplary hardware and software environment 
used to implement one or more embodiments of the inven- 
tion; and 40 

FIG. 10 schematically illustrates a typical distributed com- 
puter system using a network to connect client computers to 
server computers, according to one or more embodiments of 
the present invention. 

45 

DETAILED DESCRIPTION OF THE INVENTION 

In the following description of the preferred embodiment, 
reference is made to the accompanying drawings, which form 
a part hereof, and in which is shown by way of illustration a 50 
specific embodiment in which the invention may be practiced. 

It is to be understood that other embodiments may be utilized 
and structural changes may be made without departing from 
the scope of the present invention. 

Overview 55 

The present invention discusses the need to perform accu- 
rate dielectric and magnetic property measurements on larger 
sized samples, particularly liquids at microwave frequencies. 
These types of measurements cannot be obtained with accu- 
racy using conventional cavity perturbation methods, particu- 60 
larly for liquids, powdered, or granulated solids that require a 
surrounding container. To solve this problem, a model has 
been developed for the resonant frequency and quality factor 
of a cylindrical microwave cavity containing concentric 
cylindrical samples. This model can then be inverted to obtain 65 
the real and imaginary dielectric and magnetic constants of 
the material of interest. 
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The present invention may also determine the p of a mag- 
netic material, however, in this case, there may be typically 4 
unknowns (real and imaginary parts of both E and p), and the 
sample may be measured to isolate and determine each 
unknown separately. 

The approach is based on using exact solutions to Max- 
well’s equations for the resonant properties of a cylindrical 
microwave cavity, and also using the effective electrical con- 
ductivity of the cavity walls that is estimated from the mea- 
sured empty cavity quality factor. This new approach calcu- 
lates the complex resonant frequency and associated 
electromagnetic fields for a cylindrical microwave cavity 
with lossy walls, that is loaded with concentric, axially 
aligned, lossy dielectric and/or magnetic cylindrical samples. 
In this approach, the calculated complex resonant frequency, 
consisting of real and imaginary parts, is related to the experi- 
mentally measured quantities. 

Because this approach uses Maxwell’s equations to deter- 
mine the perturbed electromagnetic fields in the cavity with 
the material(s) inserted, one can calculate the expected wall 
losses using the fields for the loaded cavity rather titan just 
depending on the value of the fields obtained from the empty 
cavity quality factor. These additional calculations provide a 
more accurate determination of the complex dielectric con- 
stant of the material being studied. The improved approach 
may be important when working with larger samples, or 
samples with larger dielectric constants that will further per- 
turb the cavity electromagnetic fields. Also, this approach 
enables a larger sample of interest, such as a liquid or pow- 
dered or granulated solid, inside a cylindrical container. 

Technical Description 

The new technique of the present invention takes advan- 
tage of a theoretical model previously developed at JPL (Jack- 
son et al., 1 994). This more exact approach is based on exact 
solutions to Maxwell’s equations for the resonant properties 
of a cylindrical microwave cavity (Sphicopoulos et al., 1 984). 
The new feature of this measurement approach is the ability to 
use a cylindrical sample of interest, such as a liquid (like 
methane or ethane at Titan surface temperature of 92 K), 
inside a cylindrical container. This new method goes beyond 
cavity perturbation approaches (Meng et al., 1995) by allow- 
ing the diameter of a concentric cylindrical sample and holder 
to be of arbitrary sizes. 

In the mid-1 990’s, JPL had an active microwave laboratory 
where many experimental and theoretical studies were under- 
taken to understand how microwaves interact with materials. 
This approach was extended for use in earlier JPL studies on 
microwave processing of materials (Jackson et al., 1994). 
One theoretical study developed a sophisticated computer 
program to predict the temperature dependence inside a 
cylindrical sample as a function of time when a cavity reso- 
nance was excited at a given microwave power level (Jackson 
etal., 1994). 

The present invention may be used to obtain the dielectric 
permittivity and magnetic permeability of samples, by mea- 
suring the samples in a cavity with or without a sample holder. 

In one embodiment, the present invention is concerned 
with measurements made on a cylindrical sample inserted 
along the axis of a cylindrical microwave cavity using the 
TM 0n0 modes, whose resonant frequencies are independent 
of the length of the cavity, and only depend on the radius of the 
cavity as given in Eq. (1) for an ideal empty cavity. 

foM= ca O,J 2xa (!) 

where c is the speed of light in vacuum, a is the cavity radius, 
and the constant a 0 „ depends on the value n (a ol =2.405, for 
example) and is the n‘ h zero of the J 0 Bessel function. 
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The resonant frequency value of a real experimental cavity 
may be affected by any changes in the boundary conditions, 
for example due to holes in the walls or insertion of a small 
exciter or detector antenna. Of course, the resonant frequency 
may also change with temperature, since the cavity dimen- 
sions may expand or contract depending on the thermal prop- 
erties of the cavity material. 

FIG. 1(a) illustrates a simple embodiment comprising a 
cylindrical sample 100 (e.g., sample rod, or sample inside a 
rod) in a cylindrical cavity 102, situated along the cylindrical 
cavity’s axis. FIG. 1 (b) is a cross-sectional schematic of the 
sample 100 in the cavity 102 of FIG. 1(a), FIG. 1(c) is a 
cross-sectional schematic of a larger sample 100, and FIG. 
1(d) is a cross-sectional schematic of the cylindrical cavity 
102 comprising an empty cylindrical sample holder 104. 

FIG. 2(a) illustrates a cylindrical cavity 200 containing an 
annular sample 202, FIG. 2(b) is a cross-sectional schematic 
of the cylindrical cavity 200 containing an empty annulus 
holder 204, and FIG. 2(c) is a cross-sectional schematic of the 
cylindrical cavity 200 containing the sample 202 in the cylin- 
drical annulus holder 204. 

However, the present invention may be used for more gen- 
eral configurations, comprising multiple materials and 
regions. For example, in order to perform this analysis, a 
cylindrical cavity 300 may divided into concentric cylindrical 
regions 302, 304, and 306 as shown in FIG. 3 (Sphicopoulos 
et. al., 1984). 

FIG. 3 illustrates an interior of the cylindrical cavity 300 
comprising a set {i=l, 2, . . . N}(“set {i}”) of N concentric 
cylindrical zones (“the zones”) 302, 304, and 306, wherein 
the zones include a central cylinder 302 and one or more 
annular regions 304-306 surrounding the central cylinder 
302. One or more of the annular regions 304-306 may com- 
prise the sample holder, and the central cylinder 302 may 
comprise the sample or the vacuum. One of the annular zones 
(e.g., between the walls and the sample holder), e.g. 306, may 
comprise vacuum or air. A complex valued dielectric permit- 
tivity E; assigned for each zone i in the set {i} of the N zones, 
thereby forming a set of complex valued dielectric permittivi- 
ties {e ! }={e ! .=e 1 , e 2 . . . e„}. 

There may be an arbitrary number N or L of concentric 
regions. Each region 302, 304, 306 may have a different 
complex dielectric constant E 1; E 2 , E z , respectively, and 
each complex dielectric constant E 1; E 2 , E z may have a 
known temperature dependence. Given that information, the 
program (Jackson et al., 1994) may calculate the instanta- 
neous cavity complex resonant frequency at any given tem- 
perature as microwaves heat the sample. 

The program (Jackson et al., 1994) is valid for essentially 
any resonant cavity mode and any diameter sample (or con- 
centric materials) situated along the axis 308 of the cavity 
300. The outer walls 310 may be metallic walls. The outer 
walls are typically (but not necessarily) considered boundary 
conditions rather than regions 302-306. Also shown in FIG. 3 
is the radius of the cylindrical region 302, the radius p 2 of 

the cylindrical region 304, the radius p z-1 of the cylindrical 
region 306, and the radius p L of the U h or N' 7 ' cylindrical 
region. Also shown is the height or length h 312 of the cavity 
300. 

Uiis new approach of the present invention may calculate 
the complex resonant frequency and associated electromag- 
netic fields for a cylindrical microwave cavity that is loaded 
with concentric axially aligned lossy dielectric or magnetic 
cylindrical samples. The actual losses in the walls 310 are 
estimated from the measured empty cavity quality factor Q 0 . 
In this approach, the real (f ) and imaginary (f) resonance 
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frequencies f=f +jf' are related to the experimentally mea- 
sured quantities by the following expressions: 

f=f, (2) 

5 

r=f/2& o) 

where f s is the measured resonant frequency, and is the 
quality factor associated with the sample only, given by 

to 

Q,=(Wi-VQoT\ (4) 

where Q r is the loaded cavity quality factor Q. 

This calculation assumes known values for the complex 
dielectric constants of the material(s) in the cavity. However, 
15 the present invention is interested in calculating the dielectric 
constants from the measured cavity frequencies and quality 
factors, with and without the material(s) inserted. Thus, the 
present invention developed an inverse program that runs 
through a range of real and imaginary dielectric or magnetic 
20 constants until the present invention obtained the experimen- 
tal values for f s and Q^. Since this approach can actually 
determine the perturbed electromagnetic fields in the cavity 
with the material(s) inserted, the present invention can recal- 
culate the expected wall losses using the new fields rather, 
25 than just depending on the value obtained from the empty 
cavity quality factor. 

Each power sweep through the cavity resonance may be fit 
to a Lorentzian line shape with an additional linear back- 
ground frequency term. From this fit, the present invention 
30 obtained the resonant frequency, f, and quality factor, Q, for 
the measured experimental arrangement. 

The present invention has modified the program (Jackson 
et al ., 1 994) to allow the present invention to probe any range 
of real and imaginary dielectric or magnetic constants for a 
35 sample and surrounding holder, until the present invention 
obtains the complex dielectric or magnetic constants that give 
the measured cavity complex resonant frequency, f m , and 
quality factor, Q m . 

Method Steps 

40 FIG. 4 is a flowchart illustrating a method of obtaining a 
complex valued dielectric permittivity E of at least one 
sample. The method may also be used to obtain a complex 
valued magnetic permeability p of at least one sample, by 
solving for p instead of E in the discussion below. The cavity, 
45 the sample holder, and the sample may be cylindrical or 
annular regions, however, any cavity, sample, or sample 
holder that has symmetry allowing an exact solution of Max- 
well’s equations may be used. The sample may be large 
enough to perturb cavity modes of the cavity. The method 
50 may be performed for the sample in the cavity, with or without 
a sample holder. 

The method may comprise the following steps. 

Block 400 represents obtaining (e.g., measuring) one or 
more complex-valued resonance frequencies (f m ) and quality 
55 factors (Q m ) of the cavity, wherein each i m is a measurement. 
The Re(f m ) may be the measured cavity resonant frequency, 
and the Im(f m ) may be related to the measured cavity quality 
factor Q m by Im(f J=Re(fJ/2Q„,. 

The step may comprise obtaining (e.g., measuring) at least 
60 one of the f m that is a measurement of a complex -valued 
resonance frequency (f cavity ) of the cavity when empty, 
obtaining (e.g., measuring) at least one of the f m that is a 
measurement of a complex-valued resonance frequency 
(^holder) °f the cavity containing a sample holder, wherein the 
65 sample holder is without a sample or contains a vacuum, and 
obtaining (e.g., measuring) at least one of the f m that is a 
measurement of a complex-valued resonance frequency 
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(f sar „pie ) °f tlie cavity containing the sample in the sample 
holder or containing the sample without the sample holder. 

The measurements may be obtained by exciting the cavity 
with one or more cavity resonant modes. 

Block 402 represents obtaining an effective electrical con- 
ductivity ct w of the cavity’ s walls, using the i cavity . The ct w may 
be obtained using i cavity and known formulas for the calcu- 
lated complex valued f of the empty cavity (derived from 
Wilson et. al, 1949): 

°w = [Qo 2 / (^Mo/o)] (2+D ( /L) 2 /D^ (5) 

where Q 0 is empty cavity quality factor obtained from f cavity , 
p 0 is the permeability of vacuum, f 0 is the empty cavity reso- 
nant frequency also obtained from f cavity , D,. is the cavity 
diameter, and L is the cavity length. Thus, the a w is an experi- 
mentally determined electrical conductivity for the actual 
cavity walls used, because it is based on a measured i cavity , 
and can therefore take into account the real, non-ideal prop- 
erties of the cavity walls, including holes and other imperfec- 
tions in the cavity walls. 

Block 404 represents solving Maxwell’s equations exactly 
(using the method of (Sphicopoulos et. al., 1984) for a com- 
plex valued resonance frequency f of the cavity including the 
sample, wherein the f=f(E, p) is a function of the E and/or p 
of the sample. Alternatively, the Maxwell’ s equations may be 
solved for the sample in the cavity and in a sample holder, so 
that f=f(E, \i, Voider) ' s a function of both the E (and/or p of 
the sample) and the complex valued dielectric permittivity 
E h older °f the sample holder. The step may further comprise 
solving the Maxwell’s equations for a complex-valued reso- 
nance frequency g of the cavity including the sample holder 
without the sample, wherein the g=g (E. holder ) is a function of 

the E holder ■ 

More generally, the step may comprise solving the Max- 
well’s equations (Sphicopoulos et al., 1984) for a complex 
valued resonance frequency f({£„ p ; }), where a complex 
valued dielectric permittivity E, and is p, assigned for each 
zone i in a set {i=l , 2, ... N} of annular regions and a central 
cylinder, the set of zones {i} comprising the interior of the 
cylindrical cavity. The f({E„ p, }) is a single frequency for the 
whole cavity that is a function of every E t and/or p f in the 
cylinder. If there are no magnetic materials in the cavity, then 

f({e,j). 

Regions in the cavity that are air or vacuum may use the 
known complex valued permittivity for air or vacuum. For 
example, the complex permittivity for a vacuum has a real 
part equal to one and an imaginary part equal to zero. 

Block 406 represents using the ct w to modify the f to 
account for the power absorbed by the cavity’s walls, thereby 
obtaining a resultant complex valued resonant frequency F 
that is a function of the E and/or p. The step may further 
comprise modifying the g to account for the ct m „ thereby 
obtaining a resultant complex-valued resonant frequency G 
that is a function of the E holder . 

More generally, the step may further comprise modifying 
the f({E„ p,}) to account for the power absorbed by the 
cavity’s walls, thereby obtaining a resultant complex-valued 
resonant frequency F. In the general case, the modifying may 
comprise the following steps. 

(1) Calculating a surface resistance formula (e.g., p. 32, 
equation 1.125 of (Pozar, ¥ d Edition)): 

■R s =[(jtHo/te/®Ajw] 1/2 (6) 

where f r =f({E„ p,}) is the theoretical resonant frequency. 
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(2) Using R^ to obtain the power absorbed by the cavity’s 
walls, P„ a// (see e.g., (Pozar, ¥ d Edition), page 33, equation 
1.131): 

5 P^ u =PHRP-)S s \H,\ 2 ds (7) 

where t represents the tangential component of the magnetic 
field at the surface of the cavity walls, P r indicates the total 
power absorbed by the total surface area of the cavity walls, f s 
indicates an integral over the surface of the cavity walls, ds is 
to the surface element, and H ; is the magnetic field tangential 
component at the surface of the cavity walls. Thus, the energy 
loss in the cavity walls P M , a// depends on the R^ and the mag- 
netic field strength evaluated at the cavity wall for the calcu- 
lated magnetic field corresponding to the frequency f({E,., 
15 p,}). The cavity wall quality factor Q M , a// can then be deter- 
mined from P wa „ using Q wa „=2it Re(f({E,., p,))) U/P w/ 
where U is the total electromagnetic energy stored in the 
cavity (see e.g., (Pozar, ¥ d Edition), page 268, equation 6.7). 

(4) Obtaining the contribution f wa „=Re (f{E„ p,})/2 Q wa „ 
20 to the imaginary part of the resultant complex-valued reso- 
nance frequency, so that the resultant complex valued reso- 
nance frequency F({E„ p,})=f({E„ p,})+jf wall . If there are no 
magnetic materials in the cavity, then F({E,}). 

Thus, F and f may be a function of E (sample with no 
25 sample holder), or E and E holder (sample with sample holder) . 
F and f may be a function of the sample’s p if the sample is 
magnetic. 

Block 408 represents obtaining the E (or p) and E. holder by 
solving F=f /e , G=f holder , respectively, or more generally 
to by solving F({E„ \^)=i holder and F({E„ \>.,})=i sample . The 
model may choose different values for the only unknown Ej , 
to obtain a theoretical F to compare to the experimental value 
f 'sample > to see which value of E, gives a match to any desired 
degree of precision. For example, the step may comprise 
35 obtaining a range of values for the E; and selecting the value 
in the range that satisfies F=f sample to any desired degree of 
precision. The solving may be performed numerically, or 
analytically, for example. 

The cavity may comprise one or more additional materials 
40 (e.g., sample holder, vacuum, air, etc.) having known material 
dielectric permittivity and known material magnetic perme- 
ability, such that the E is an only unknown during the solving 
step. In this regard, the method may first obtain the dielectric 
properties of the unknown materials in the cavity. For 
45 example, the present invention may obtain the Voider by 
solving G=f hoi der \ and then solve of F=f sample for the sample 
inserted in the sample holder, using the G: holder obtained by 
the solving of G=i hoIder . 

Thus, more generally, the step may comprise obtaining the 
50 complex-valued dielectric permittivity of an empty sample 
holder E holder in a zone i=I, or in two non-adjacent zones I and 
K occupied by the holder, where the dielectric permittivity E t 
(and p, if necessary) for each of the other zones in the set {i} 
is known, by solving F({E ( , p,})=f m , where f m is the complex- 
55 valued resonance frequency for the cavity containing an 
empty sample holder. 

The step may then further comprise obtaining the complex- 
valued dielectric permittivity E or p of a sample, in a zone i=J, 
contained in a sample holder, where the complex -valued per- 
60 mittivity E,(and p, if necessary) for each of the other zones in 
the set {i} is known, by solving F({E ; ., p,} )=f,„, where f m is the 
complex-valued resonance frequency for the cavity contain- 
ing a sample in a sample holder. The E hoIder solved in the 
previous paragraph may be used as one of the known E t . 

65 Thus, Blocks 404-408 above represent solving Maxwell’s 
equations exactly for E and/or p and/or E holder using the f m as 
known quantities, thereby obtaining the E and/or p of the 



US 8,653,819 B2 


11 

sample and/or ^ holder of the sample holder. The method of the 
present invention may comprise exciting the cavity, measur- 
ing quality factors and resonant frequencies to obtain the 
experimental complex resonant frequency f of the cavity, 
solving Maxwell’s equations for the theoretical complex 5 
resonant frequency F, and obtaining E and/or p by equating 
the theoretical complex resonant frequency F to the measured 
complex resonant frequency f. 

The present invention allows for exploring any range of the 
unknown complex E or to determine the value that satisfies to 
F=f to any degree of precision desired. 

Block 410 represents using the E or p to interpret remote 
sensing data obtained from earth or extra-terrestrial material, 
thereby identifying one or more compositions, optical prop- 
erties, or dimensions of the material; and transmitting, stor- 15 
ing, or displaying the compositions, optical properties, or 
dimensions of the material (Janssen et al., 2008). This inven- 
tion could also support data interpretation from future flight 
missions including SMAP (Soil Moisture Active & Passive) 
(1.2-1. 4 GHz), and DESDynl (Deformation, Ecosystem 20 
Structure and Dynamics of Ice) (1.2 GHz). 

However, the present invention is not limited to providing 
E or p for any one particular application. The method of the 
present invention can be used to provides or p for any appli- 
cation. 25 

In general, the present invention solves for one unknown at 
a time. For the case of non-magnetic materials in the cavity, 
the functions f and F are not a function of p In the case of 
magnetic materials in the cavity, the permeability of the mate- 
rials in the cavity will have to be known. The sample’s per- 30 
meability p can be determined using the method of the present 
invention once the sample’s permittivity is known, as further 
discussed below. 

Obtaining Magnetic Permeability 

A similar method may be used to measure magnetic per- 35 
meability, as illustrated in FIG. 5. The present description 
given above may make accurate permittivity measurements 
on essentially any size non-magnetic materials. 

However, even if the sample is magnetic, the permittivity 
may be determined by placing a small cylindrical sample and 40 
holder along the axis of a cylindrical cavity where the mag- 
netic field is close to zero. 

For example, Block 500 of FIG. 5 represents determining 
the E of a magnetic sample. The step may comprise using a 
cavity resonant mode that is at least one Transverse Magnetic 45 
(TM 0 „ 0 ) mode. The step may further comprise performing the 
steps of Blocks 400-408 in FIG. 4 for the sample that is a first 
cylindrical sample positioned along the cavity’s axis, the first 
cylindrical sample having a first radius 106 that is sufficiently 
small so that the magnetic field’s strength (produced by the 50 
cavity resonant mode) at the first radius’ end 108 is negligible 
(see also FIG. 1(a)), thereby allowing the Maxwell’s equa- 
tions and the F to be solved for only one unknown, the E. 

The present invention may also measure the permeability 
of samples. The permeability can be determined by first using 55 
a thin sample to obtain the complex dielectric constant as 
mentioned above (e.g., using the first cylindrical sample). 
Then a larger diameter sample (second cylindrical sample) is 
used so it sees both a magnetic and electric field produced by 
the cavity resonant mode. Knowing the dielectric constant, 60 
the present invention can then back out the permeability. 

Block 502 represents determining the p of a magnetic 
sample. The step may comprise performing the steps of 
Blocks 400-408 in FIG. 4, using the sample that is the second 
cylinder comprising a same material as the first cylinder, and 65 
having a larger radius than the first radius, such that the 
second cylinder positioned along the axis experiences both 
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the magnetic field and the electric field, thereby allowing the 
Maxwell’s equations and the F to be solved for only one 
unknown, the p, using the E obtained using Block 410. 

With this approach there may be a slight shift in the reso- 
nant frequency between the small and large samples. With a 
proper design of the cavity, these permittivity and permeabil- 
ity measurements may be made for at least the three lowest 
resonant frequencies (TM 010 , TM 020 , and TM 030 ). 

The present invention further discloses a second method 
for measuring the permeability of a sample, where only one 
sample is required. In this case, the cavity resonant modes 
may comprise a first frequency, a first magnetic field, and a 
first electric field associated with a first TM 0 „ 0 mode, where n 
is an integer. The cavity resonant modes may further comprise 
a second frequency, a second magnetic field and a second 
electric field associated with a second TM 0 „ 0 mode, wherein 
the first frequency and the second frequency such that the E 
and the p at the first frequency are negligibly different from 
the E and the p at the second frequency. 

In this case, Block 500 represents (1) exciting the cavity 
with the first TM 0 „ 0 mode where the sample is an annulus 
positioned, placed, or located at a radius value of the cavity 
associated with a maximum of the first electric field and a zero 
of the first magnetic field (where the magnetic field is zero); 
and (2) performing the steps of Blocks 400-408 in FIG. 4 
using the first TM 0 „ 0 mode, thereby allowing the Maxwell’s 
equations and the F to be solved for only one unknown, the E. 

Block 502 then represents performing the steps of Blocks 
400-408 in FIG. 4 using the second TM 0 „ 0 mode so that the 
sample experiences both the second electric field and the 
second magnetic field, thereby allowing the Maxwell’s equa- 
tions and the F to be solved for only one unknown, the p, using 
the E obtained using the first TM 0n0 mode. 

In one example, the first TM mode is a TM 020 mode, 
measurements for this TM 020 mode may determine the per- 
mittivity, and the second TM mode is a TM 010 mode that 
allows the same sample to be in both magnetic and electric 
fields, thereby allowing the present invention to back out the 
permeability knowing the permittivity obtained using the 
TM 020 mode. The annular sample may a solid, or a liquid or 
powder in an annular holder, for example. 

The present invention has developed Mathematical pro- 
grams for calculating the permittivity and permeability using 
this technique. These programs are able to determine the 
changes in the electromagnetic field at the cavity wall from its 
empty cavity and loaded cavity states. This ability goes 
beyond the previous perturbation approaches. 

Block 504 represents using the p, E. 

Rootfinder Method For Solving F({E ; ., p,})-f„=0 

A highly efficient rootfinder that may be used to solve an 
equation such as F({E„ p ; })-f m =0 for a complex valued 
dielectric permittivity E z or complex valued dielectric perme- 
ability p, has been developed and is described below. This 
rootfinder may be applied to a large variety of problems. 

One or more solutions of the equation C(z)=0, where C=Re 
C+j Im C and z=Re z+j Im z, can be obtained as intersection 
points of two lines that lie in the complex z-plane. The coor- 
dinate axes in that plane are labeled Re z and Im z. The two 
null, or zero, lines that intersect can be determined by the 
conditions Re C(z)=0, Im C(z)=0. The procedure for locating 
the one or more solutions of C(z)=0 may be made highly 
efficient by first obtaining approximate solutions with a mini- 
mum of calculation. Only then are complicated quantities 
evaluated that accurately locate the solution or solutions. 

FIG. 6 is a flowchart illustrating the rootfinder method. 
FIG. 7 illustrates a mesh in the complex plane used in the 
method. 
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The first step in the procedure, represented in Block 600, is 
to specify a range of values of z, within a rectangle 700 in the 
complex z-plane 702, that is to be searched for solutions of 
C(z)=0. The rectangle 700 is specified to be the region having 
coordinates along the Re z axis between ZRMIN and 5 
ZRMAX, and coordinates along the Im z axis between 
ZIMIN and ZIMAX. 

Next a mesh is introduced in this rectangle, as illustrated in 
Block 602, the mesh having NR equal intervals 704 between 
ZRMIN and ZRMAX, and NI equal intervals between to 
ZIMIN and ZIMAX. Each small rectangle 706 in the mesh 
will be associated with the indices (L, M) that locate its lower 
left-hand comer. 

Block 604 represents scanning and evaluating Re C(z) on 
the mesh. Start at L=l, M=l, scan from left to right, and then 15 
move upward to the next row and scan from left to right. 
Repeat this procedure until the entire large rectangle has been 
scanned. Evaluate the function Re C(z) at each comer of each 
rectangle 706 in the mesh, i.e. at each mesh point 708. Then 
return to L=l, M=1 and scan again in the same manner, but 20 
also determine sides 710 of the small rectangles 706 in the 
mesh such that Re C(z) changes sign there, as represented in 
Block 606. For each rectangle (L, M) assign the value 0 or 2 
to a function JR(L,M) according to whether there are zero or 
two sides 710 of the rectangle 706 where Re C(z) changes 25 
sign. 

A similar procedure is carried out for Im C(z)=0, as illus- 
trated in Blocks 608 and 610, but now assign the value 0 or 2 
to a function JI(L, M) according to whether there are zero or 
two sides 712 of the small rectangle 706 where Im C(z) 30 
changes sign. 

The next step, represented in Block 612, is to scan the large 
rectangle 700 and select for further processing each small 
rectangle 714 (thick dashed rectangle in FIG. 7) where JR(L, 
M)=2 and JI(L, M)=2. Any point z within or on the boundary 35 
of such a rectangle (L, M) 714 is an approximate solution of 
C(z)=0. However a more accurate solution can be obtained by 
the following steps. 

In each rectangle (L, M) 714 that satisfies the conditions 
specified, points 716a, 7166 on the sides 710, 712 of the 40 
rectangle 714 such that Re C(z)=0 can be obtained to good 
approximation by linear interpolation of values of Re C(z) on 
the comers of the rectangle 714, as represented in Block 614. 
The two null points 716a, 7166 on sides of the rectangle (L, 

M) 714 define a straight null line segment for the fiinction Re 45 
C(z). This line segment lies in the rectangle (L, M) 714. 

For that same (L, M) rectangle 714, points 718a, 7186 on 
the sides of the rectangle such that to good approximation Im 
C(z)=0 can be obtained by linear interpolation of Im G(z) on 
the corners of the rectangle (Block 616). Those two null 50 
points 718a, 7186 define a straight null line segment for the 
function Im C(z). That line segment lies in the rectangle (L, 

M) 714. 

The intersection point z 0 of the two straight null lines inside 
the rectangle (L, M) 714 can be found by simple algebra, as 55 
represented in Block 618. The value of z 0 is a better approxi- 
mate solution of C(z)=0. 

The entire procedure of Blocks 600-618 is then repeated, as 
represented in Block 620, by further partitioning the above 
rectangle (L, M) 714 with a smaller mesh size and determin- 60 
ing a new intersection point, which will now be called z r 

Hie difference [Re (Zj)-Re (z 0 )] and the difference [Im 
(zj-lm (z 0 )] are now individually compared with a specified 
precision parameter, as represented in Block 622. If each of 
these differences meets the specified precision criterion, then 65 
Zj will be taken as the solution of C(z)=0. If the precision 
criterion is not met, then this procedure is repeated with 
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successively smaller rectangles covered by successively finer 
meshes until the specified precision criterion is met by two 
successive calculations. 

Block 624 represents the final value for the intersection 
point, now called z^, is a solution of C(z)=0 that meets the 
precision that has been specified. Application of this proce- 
dure an infinite number of times will converge to an exact 
solution of C(z)=0. 

Thus, as applied to the solving of F=f sample , the method of 
the present invention may further comprise finding roots of 

F -U mp / e =0. U: 

(i) specifying a range of values for z within a first rectangle 
in a complex plane, wherein the z are complex numbers in a 
complex plane (Block 600); 

(ii) introducing a mesh of second rectangles into the first 
rectangle, thereby forming the mesh comprising mesh points 
at corners of second rectangles, wherein the second rect- 
angles are smaller than the first rectangles (Block 602); 

(iii) scanning and evaluating Re F (z) on the mesh at each of 
the mesh points (Block 604); 

(iv) determining one or more first sides of the second 
rectangles where the Re (F(z)) changes sign (Block 606); 

(v) scanning evaluating Im F(z) at each of the mesh points 
(Block 608); 

(vi) determining one or more second sides of the second 
rectangles where the Im F(£) changes sign (Block 610); 

(vii) selecting the second rectangles having both the first 
sides and second sides, where both the Re F (z) and the Im F (z) 
change sign, to obtain selected second rectangles (Block 
612); 

(viii) interpolating the Re F(z) in the selected second rect- 
angles to define a Re F(z) null line (Re C(z)=Re F(z)=0) 
(Block 614); 

(ix) interpolating the Im F(z) in the selected second rect- 
angles to define an Im F(z) null line (Im C(z)=Im F(z)=0) 
(Block 616); 

(x) finding a first intersection point between the Re F(z) 
null line and the Im F(z) null line (Block 618); 

(xi) setting the second rectangle as the first rectangle and 
repeating steps (i)-(x) using the second rectangle as the first 
rectangle, thereby finding a new or second intersection point 
(Block 622); 

(xii) obtaining a difference between the first intersection 
point and the second intersection point, and comparing the 
difference with a precision parameter that is equal to a desired 
precision (Block 622); 

(xiii) repeating the method until the difference is at least as 
small as the precision parameter, thereby obtaining a final 
intersection point; 

(xiv) setting the final intersection point equal to the prop- 
erty (Block 624). 

Steps may be omitted or added as desired. 

Experimental Results 

Polyvinyl Chloride (PVC) Rod Samples 

In one embodiment, the present invention used an existing 
microwave cavity (3" long cylindrical cavity) designed to 
excite the TM 010 resonant mode with a resonant frequency of 
2.45 GHz (using equation (1), the TM 0 „ 0 mode has 
f 010 =2.452 GHz). The test samples were Polyvinyl Chloride 
(PVC) rods, some of which were enclosed in a Polycarbonate 
(PC) tube. Hie present invention performed dielectric prop- 
erty measurements on various sized PVC rods (Type 1 PVC 
solid rods of Va” diameter and %" diameter, simulating e.g., 
the sample), either alone or enclosed in the PC tube (used, e.g. 
as the sample holder). A 3 /s" outer diameter (OD)x!/i6" wall 
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PC tube was used as the holder for the smaller 14" diameter 
PVC rod (that could also simulate a liquid contained in the 
holder). 

These measurements should give the same values for the 
PVC real and imaginary dielectric constants for both PVC 5 
rods. 

The present invention was careful regarding the way 
energy is transmitted into and out of the cavity. The cavity was 
excited using a small loop antenna positioned to couple to the 
TM 010 mode magnetic field. The energy in the cavity was to 
detected using a similar small loop antenna situated on the 
sidewall. The approach taken by the present invention was to 
weakly couple the microwave energy into the cavity. This was 
accomplished by designing a special support for a rotatable 
rod attached to the loop antenna. The rod was rotated so that 15 
the excitation loop area barely coupled to the magnetic field 
lines associated with the TM 010 mode. This approach requires 
a much larger drive signal, since -99% of the signal may be 
reflected back. In this way, a negligible amount of the energy 
that does enter the cavity may leak back out through the 20 
weakly coupled exciter or detector loops, and cause an error 
in the quality factor (Q) measurement. 

An Agilent Technologies E8364B Vector Network Ana- 
lyzer (VNA) was used to sweep the frequency through the 
resonance. With the exceptional signal to noise characteris- 25 
tics of this unit, only 8 mW of drive power was needed to 
obtain measurable signals. 

First the empty cavity resonant frequency and quality fac- 
tor were measured. These empty cavity measurements were 
needed because the cavity walls absorb microwave energy in 30 
addition to the sample and a correction for this additional 
energy loss needed to be taken into account. 

Then, the sample of interest was inserted, and the loaded 
sample cavity frequency and quality factor were measured. 

For the case of the 14" diameter PVC rod inserted into the PC 35 
tube, the PC tube complex dielectric constant was first deter- 
mined by inserting this tube into the cavity separately, and 
making the frequency and quality factor measurements for 
the tube alone. Then, knowing the dielectric properties of the 
tube, the present invention repeated the measurements with 40 
the PVC rod inserted in the PC tube. The present invention 
also separately performed cavity measurements on the larger 
Vs" OD PVC rod. 

Each power sweep through the cavity resonance may be fit 
to a Lorentzian line shape with an additional linear back- 45 
ground frequency term. From the fit, the present invention 
obtained the resonant frequency, f o =2.38071e+09+/-12 Hz, 
and quality factor, Q 0 =1414.27+/-0.03, for one example of a 
measured experimental arrangement. 

These sets of measurements were used to determine the 50 
dielectric properties of the two different sized PVC rods using 
various calculation methods (cavity perturbation and the 
present invention). The validation of the new technique of the 
present invention depended on the measured values of the two 
PVC rods being consistent with each other. 55 

The difference in the real and imaginary dielectric con- 
stants between the two PVC rods was 3.6% and 1.1%, respec- 
tively, using the approach of the present invention. These 
results confirmed the validity of the present invention’s new 
technique. The resultant deviations were partially due to the 60 
slight non-uniformity in the PVC rods, the fact that there was 
a slight gap between the 14" PVC rod inserted in the PC tube. 

Ice Samples 

This technique was also validated for a solid ice sample 
enclosed in a container at lower temperatures. 65 

Additional measurements were made on pure water ice 
samples at -20° C. and -38° C. For the real part of the 
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dielectric constant at these temperatures, the present inven- 
tion obtained 2.91 and 2.99 respectively. The only other reli- 
able measured values of the real part of the dielectric constant 
in the present invention’s frequency and temperature ranges 
(0.15 to 2.5 GHz, -1 to -60° C.) have been obtained by 
Westphal (Ulaby), and are in the range 2.90-2.95. These val- 
ues by Westphal are in excellent agreement with the present 
invention’s measurements at -20.3° C. and -38.0° C. 

However, a comparison of the present invention’s imagi- 
nary dielectric constant measurements of pure water ice was 
difficult because there is substantial scatter in the previously 
published studies (Ulaby et. al, 1986). Westphal (c.f. Warren 
et al ., 1 985) noticed that annealing ice samples for a few hours 
at - 10° C. resulted in lower losses, and used this technique for 
making the ice measurements. However, even with this 
approach, Westphal estimated the accuracy in measuring the 
imaginary dielectric constant to be ~+20%. 

The present invention always cooled the cavity as quickly 
as possible to the target low temperatures. However, because 
of Westphal’s uncertainties, the present invention was not 
able to make a meaningful comparison of the present inven- 
tion’s imaginary dielectric constant with literature values. 
Nevertheless, the results of the present invention’s initial 
study of PVC and ice samples gives confidence in using the 
present invention’s new approach for performing further 
cryogenic measurements on liquid methane and ethane. 

Condensed Liquefied And Cryogenic Samples 

One of the motivations for developing this invention was to 
be able to perform accurate dielectric measurements on cryo- 
genic liquids, such as methane and ethane. These are materi- 
als assumed to be components in lakes on Titan, a moon of 
Saturn. The Cassini spacecraft has been taking radar mea- 
surements of Titan’s lake at a frequency of 13.78 GHz. In the 
past, the “free-space,” “resonant-cavity,” “capacitive,” and 
“calorimetric” methods (Whiffen, 1950; Lorenz, 1999; 
Lorenz and Shandera, 2001 ; Lorenz et al., 2003; Paillou et al., 
2008; Paillou et al., 2008a) have been used to measure the 
complex dielectric constant of some Titan-relevant materials. 
Some of the early studies were performed at lower frequen- 
cies, near room temperature and were not very accurate 
(Whiffen, 1950; Lorenz and Shandera, 2001; Lorenz et al., 
2003). More recent studies by Paillou and co-workers (Pail- 
lou et al., 2008; Paillou et al., 2008a), working in the micro- 
wave region, have used the free-space method (Ghodgaonkar, 
et al., 1989) to measure the complex dielectric constant of 
cryogenic liquids and the resonant-cavity method to measure 
solids. The more accurate of these earlier resonant-cavity 
studies have used a form of the cavity perturbation to extract 
the dielectric constants from their microwave measurements 
(Vaid et al., 1 979; Rodriguez et al., 2003). The present inven- 
tion has fabricated a cylindrical microwave cavity designed to 
excite the TM 0 „ 0 modes in order to perform these cryogenic 
microwave measurements. By making the cavity length 
small, the higher order unwanted modes are moved to higher 
frequencies, allowing the first three TMo„ 0 modes to be unper- 
turbed. 

The radius for an empty cavity tuned to have a TM 020 
resonant frequency at 13.78 GHz is a -1.91 cm. However, 
when a quartz tube used as a sample holder is inserted in the 
cavity, the resonant frequency may be lowered. For the reso- 
nant frequency to be at 13.78 GHz when a thin-walled quartz 
tube filled with liquid methane is inserted, the present inven- 
tion estimates the cavity radius to be ~1 .65 cm. 

A drawing of the custom microwave resonant microwave 
cavity 800 fabricated for measuring condensed liquids is 
shown in FIG. 8. The cavity is sealed to eliminate water 
condensation. The cavity may be first evacuated and then 
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filled with a small amount of transfer gas to provide good 
thermal contact between the quartz tube 802 and cavity ’ s 800 
copper walls. The quartz tube 802 may be moved in and out of 
the cavity 800 and may be sealed by a Teflon O-ring. A gas 
handling system may be used to allow pure hydrocarbons or 5 
their mixtures (such as methane and ethane) to be prepared 
and condensed through a capillary fill line 804 into the quartz 
tube 802 at cryogenic temperatures. A pressure gauge may 
measure the pressure in the transfer line. Once liquid forms, 
the measured pressure corresponds to the vapor pressure, to 
which is a function of the interface temperature. The experi- 
menter knows when the quartz tube 802 is filled, since at that 
time, the gas-liquid interface may move up the wanner fill 
line 804 , where a higher vapor pressure may be measured. 
Also shown is the exciter 806 (wherein exciting radiation is 15 
inputted) and the detector 408 (for detecting the resonance 
frequency and Q factor, e.g., with a loop antenna) 

A larger microwave cavity having the same design as 
shown in FIG. 8 may be fabricated to operate at 2.38 GHz. 
This new cavity may be used to study solid icy compositions 20 
relevant to understanding the surface behavior of Saturn’s 
airless icy satellites. A different quartz tube design may be 
needed when solid samples are to be initially inserted into the 
quartz tube and measured. These solid samples may be char- 
acterized, before and after their dielectric properties are mea- 25 
sured, using a cryogenic polarizing microscope, and low tem- 
perature scanning electron microscope with an X-ray Energy 
Dispersive Spectrometer (EDS). 

A cryogenic facility may be constructed to provide preci- 
sion temperature control of the disk-shaped microwave cav- 30 
ity. The cavity may be supported inside a dewar using thin 
walled stainless steel (SS) tubes. The dewar may be filled with 
liquid nitrogen to a level below the cavity. A long copper rod 
or finger may extend to near the bottom of the dewar to act as 
a cooling link to the cavity. A stainless steel insert, inserted 35 
between the copper rod and cavity may act as a thermal 
impedance. With this design, the cooling power to the cavity 
should remain rather constant, even as the liquid nitrogen 
level is reduced. The cavity may have a heater and thermom- 
eter, and the cavity temperature may be controlled using a 40 
proportional, integral, differential (PID) controller. The tem- 
perature of the cavity may have a stability better than 0.01 ° C. 
The cavity can be easily detached as a module if there is a 
need to insert a solid sample into the quartz tube. The tem- 
perature controllers and vacuum pumps needed for this cryo- 45 
genic system may also be provided, as well as a VNA (such as 
an HP 85 10C VNA) for measuring the resonant shape. 

There are three steps in the procedure for performing mea- 
surements. First, the resonant frequency and quality factor of 
the empty cavity may be measured over the temperature range 50 
of interest. This is needed to determine the temperature 
dependence of the cavity wall electrical conductivity, from 
which the present invention can obtain the absorption behav- 
ior of the cavity walls. Then, the empty quartz tube may be 
inserted into the cavity and the measurements repeated. With 55 
these two sets of initial measurements, the present invention 
may calculate the complex dielectric constant of the quartz 
tube over the measured temperature range. Then, the final set 
of measurements may be made with the sample and quartz 
tube inserted. With these data sets, the present invention may 60 
calculate the complex dielectric constant of the sample over 
the entire temperature range of interest. 

Error Analysis 

The absolute accuracy in the calculated dielectric constants 
primarily depends on the accuracy in (1 ) the measurement of 65 
the microwave cavity and quartz tube dimensions, (2) the 
measurement of the shift in the cavity resonant frequency and 
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quality factor when the sample is inserted, and (3) the calcu- 
lation method. The microwave cavity and quartz tube dimen- 
sions, as well as the frequency shift, can be measured quite 
accurately. The new calculation method depends on an exact 
solution of Maxwell’s equations for this experimental con- 
figuration. However, it assumes perfect boundary conditions. 
In reality, there can be perturbations in the electromagnetic 
fields due to no ideal boundary conditions (holes in cavity 
walls, exciter and detector antennas and etc.) that modify the 
currents in the cavity walls. These perturbations may be mini- 
mized by using a metal coating on the top of the inserted 
quartz tube that may approximate the expected metal bound- 
ary, and may also significantly reduce the leakage of energy 
out the cavity opening. The present invention may also use 
specially designed antennas with very small metal loops that 
approximate the expected boundary. 

With regards to determining the complex dielectric con- 
stant (where E' is Re (£), E" is Im(E), and the E" may be 
determined from the measured quality factor shift), there may 
be a limitation as to how small a E" can be determined. For the 
worst case scenario of a pure methane sample with a small 
loss tangent of, say tan 5=E7E'=1CT 5 , the present invention 
estimates the shift in the cavity Q, when the sample is 
inserted, may only be -0.35%. For an expected initial Q of 
~1 1 ,300, this corresponds to a shift in Q of ~40. This accuracy 
may depend on having a good signal to noise sweep measure- 
ment, and on how well the measured resonant line shape 
matches the expected Lorentzian shape. The present inven- 
tion may be able to measure this low loss tangent to ~5%. 

Possible Modifications 

The f„, and Q„, may be measured by exciting the cylindrical 
cavity with a TM 0 „ 0 mode of microwave radiation, with n 
such that (but not limited to) 0<n<3. However, the present 
invention is not limited to a particular excitation frequency. 

The cavity may have any size, for example a radius of the 
cavity may be selected to excite the TM modes for the 
required frequency range. 

The f„, and the Q m may be measured or obtained as a 
function of temperature, thereby obtaining E, E. holder as a 
function of temperature. 

Advantages and Improvements 

The present invention allows accurate measurements of 
dielectric permittivity and magnetic permeability of various 
materials. The real part dielectric permittivity is a measure of 
the reflectivity of the material, for example, and the imaginary 
part of the dielectric permittivity is a measure of the absorp- 
tion of the material, for example. 

The present invention uses one or more modes or reso- 
nances of a cavity, for the empty cavity and for the cavity 
containing a sample and/or sample holder. The sample and 
the sample holder perturb the empty cavity modes and shift 
the field behavior. The shift depends on the size of the sample. 
Prior to the present invention, only permittivities for samples 
having small sizes of the order of a few percent of the cavity ’ s 
volume could be measured accurately, using a cavity pertur- 
bation approach. 

Prior to the present invention, attempts to measure the 
change in Q factor with and without the sample measured 
only small samples (a few percent of cavity volume), ignored 
certain errors, and ignored the effect of the sample holder, and 
used only solid samples. 

The present invention, on the other hand, may use a larger 
sample, and accounts for cases where the fields at the cavity 
walls are changed due to the presence of the sample, as 
compared to the empty cavity. Therefore, the present inven- 
tion may determine the new field at the walls when the sample 
gets larger. The present invention may also compensate for 
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imperfections (losses) in the cavity walls. The present inven- 
tion may use a rod of solid or a liquid in a sample holder. 

The present invention may solve Maxwell’s equations 
exactly for a cylindrical cavity, containing a sample that is a 
cylindrical tube, or containing the sample that is a cylindrical 5 
tube within the sample holder that is also a cylindrical tube. 
However, the present invention’s cavity may be a cylinder of 
any size that can be broken up into an arbitrary number of 
regions or cylinders, e.g., 20, and solve for the behavior in 
each region. The present invention may monitor the tempera- to 
ture profile inside the sample rod as a function of time, for 
different materials, thereby obtaining the properties of the 
sample and sample holder as a function of time and tempera- 
ture. For example, the present invention may account for the 
inside of the sample heating before the surface of the sample. 1 5 
Using the present invention, the permittivity of the sample 
obtained with and without the sample holder is the same to 
within a few percent. Therefore, the present invention can 
more accurately account for the effect of the sample holder on 
the fields. 20 

The validation of the present invention’s new approach 
opens up opportunities for performing complex dielectric 
constant measurements on larger sized samples (both solids 
and liquids and porous materials) that camiot accurately be 
obtained by other means. An important application of this 25 
new approach is the determination of the dielectric properties 
of candidate liquid hydrocarbon compositions (methane/ 
ethane) that may be the components of the apparent lakes 
observed on Titan (~90 K). These new measurements may 
support the interpretation of radar echoes obtained from the 30 
Cassini Synthetic Aperture Radar Imager (SAR) operating at 
13.8 GHz. 

These types of low temperature measurements may require 
new cavity designs and a better temperature controlled cryo- 
genic chamber. Small disk shaped cylindrical cavities can be 35 
designed to excite the first three TM 0 „ 0 modes without inter- 
ference from other cavity modes, i.e., foio 2 GHz, f 020 5 
GHz, and f 03 o = 8 GHz. The frequency range from -1-15 GHz 
could be covered by using two different sized chambers. 

Examples Of Applications 40 

Applications of the present invention described herein are 
merely provided as examples. 

The present invention may perform cryogenic microwave 
dielectric property measurements to support the interpreta- 
tion of radar and radiometric observations of Titan and other 45 
icy satellites of Saturn. In particular, the present invention 
may use these measurements to interpret Cassini radiometric 
and radar observations at 13.78 GHz (2.2 cm) of Titan, and 
other Saturnian icy satellites obtained by the Cassini RADAR 
instrument, and available from the Cassini Planetary Data 50 
System (PDS) archive, as well as ground-based radar mea- 
surements collected on the main airless icy satellites of Saturn 
at 2.38 GHz (13 cm). New dielectric measurements are essen- 
tial to constrain the surface and near-surface composition of 
these bodies. 55 

In one embodiment, there are four main tasks. 

(1 ) Setting up a microwave cavity and cryogenic system for 
operation at temperatures down to near liquid nitrogen (-80 
K). This first microwave cavity may make measurements of 
complex dielectric properties at -13.78 GHz (2.2 cm) and 60 
have the capability of measuring low loss tangents to a high 
degree of precision (e.g. down to -1 0 -5 to within -5%) using 
the new more accurate analysis model. 

(2) Subsequently measuring properties of liquid hydrocar- 
bons relevant to Titan, specifically methane, ethane and their 65 
mixtures initially without, and afterwards with impurities 
(dissolved nitrogen, solid hydrocarbons, etc.). 
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(3) Applying these results to the analysis of lakes on Titan. 
The present invention may estimate lake depth observation 
limits as a function of composition, and RADAR looks and 
geometry. These dielectric measurements may be used in 
conjunction with RADAR stereo bathymetric measurements 
of apparent subsurface intra-lake structures, to (a) test 
whether these structures are consistent with being observed 
through volumes of cryogenic liquids, (b) bracket compatible 
composition ranges using a wide range of dielectric proper- 
ties measured in the lab, and (c) refine and apply a technique 
to estimate the depth of subsurface intra-lake structures if 
viable dielectric properties are in a narrow range. 

(4) Performing additional dielectric measurements on sol- 
ids, mainly various ice mixtures, to support the interpretation 
of the publicly available microwave radar and radiometry 
observations of the main Saturnian icy satellites, namely 
Titan, Mimas, Phoebe, Enceladus, Tethys, Dione, Rhea, 
Hyperion and Iapetus. For some of these studies, different 
microwave cavities operating at various frequencies may be 
fabricated. 

Expected Significance 

Major objectives of planetary research revolve around 
understanding the origin and evolution of solar system 
objects. In particular, the surfaces of the outer icy satellites 
contain a record of their history. Constraining their composi- 
tion is necessary to determine the physical or geological 
endogenic or exogenic processes that have affected their sur- 
face, and to determine the state and early evolution of the solar 
system, since the materials present at the surface of these 
satellites are those from which the solar system formed. 

Radar and radiometry observations can make an important 
contribution to this effort. However, the interpretation of the 
accumulating data requires knowledge of the properties of the 
materials that comprise these surfaces. Existing data is sparse 
on properties of candidate surface materials at the low tem- 
peratures on these satellites, particularly dielectric properties 
at microwave frequencies. The first part of the proposed effort 
will fill gaps and close uncertainties in these properties of one 
of the most important of these, namely, liquid methane- 
ethane mixtures, thereby enabling more thorough interpreta- 
tion of data obtained by the Cassini RADAR instrument on 
the lakes observed on Titan. Additional dielectric property 
measurements may be performed on ice mixed with various 
contaminants (annnonia, silicate, organics, etc.) and organic 
compounds, at temperatures relevant to the main icy satellites 
ofSaturn(~60-100 K) and at other frequencies. Most of these 
measurements will be totally new and are critical to better 
understand the surface composition of Saturn's icy worlds. 

Lines of Inquiry 

I. Titan’s Hydrocarbon Lakes and Seas 

a) Observation of Titan’s Lakes and Seas 

Since 2004, Titan has been explored extensively at 2.2 cm 
(13.8 GHz) wavelength by the RADAR instrument through 
the Cassini prime mission, in scatterometery and altimetry 
modes, as well in the high-resolution synthetic aperture radar 
(SAR) mode (Elachi et al., 2005; Lunine et al. 2008; Wye et 
al., 2007). Radiometric data at the same wavelength were 
obtained concurrently in all radar modes (Janssen etal. 2009). 

The observations have been interpreted to yield unique 
information about Titan’s diverse and complex surface, 
which includes vast dune fields, and evidence for the shaping 
and evolution of the surface by fluvial erosion, putative cryo- 
volcanism, impacts, and tectonics (Lorenz et al., 2006; Stofan 
et al., 2007; Lorenz et al., 2008; Le Mouelic et al., 2008; Wall 
et al., 2009). In particular, strong evidence has been found for 
liquid-filled seas and lakes on Titan, especially in the high 
northern latitudes (Stofan et al., 2007). The recent detection 
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of liquid ethane in Ontario Lacus by the VIMS instrument 
(Visual and Infrared Mapping Spectrometer) (Barnes et al., 
2009) confirmed a long-held idea that these lakes and seas are 
presently filled with liquid (Stofan et al., 2007; Hayes et al., 
2008; Mitchell et al., in review), most likely primarily a 5 
combination of methane (due to a methane-rich atmosphere, 
VIMS may not detect surface methane) and ethane. 

An accurate measurement of the microwave loss tangent of 
liquid ethane and methane and their mixture is of prime 
interest since it will allow the present invention to obtain a 10 
reliable estimate of the depth of the lakes. For example, the 
-2x10 “ 3 loss tangent obtained by Paillou et al. (2008) based 
on a Liquefied Natural Gas (LNG) mixture (primarily con- 
taining methane, but also ethane, and various other volatiles 15 
and hydrocarbons, including nitrogen), implies rather shal- 
low (up to -7 m) lakes where the bottom is even slightly 
visible. This is not inconsistent with observations of lake 
Ontario Lacus, which exhibits SAR backscatter that falls off 
exponentially with distance from the margins. When one 20 
extrapolates local altimeter data with slopes of ~ 10 -3 , the 
observations would appear to suggest a loss tangent of ~7.6x 
10 -4 (Hayes et al., in prep.), differing from Paillou et al.’s 
(2008) measurement by a factor of less than two. 

However, this contrasts with more direct stereo bathymet- 25 
ric measurements of the depth of seas in the Arctic calculated 
from stereo mappings (Digital Topographic Mapping, or 
DTM) of the polar regions led by collaborator R. Kirk (Kirk 
et al., 2010). A surprising result has been the apparent success 
of the automatic image-matching algorithm in very dark por- 11 1 
tions of the seas, though not in the darkest and completely 
featureless areas. The DTM indicates a few hundred meters of 
relief in areas believed to be liquid-covered. Some of the 
DTM results in the darker parts of the seas may well result , 
from spurious matching between speckle noise in the images, 
but substantial portions with backscatter cross-sections 
between -10 and -15 db (e.g., around the northern sides of 
Mayda Insula) clearly show channels and other features com- 
mon to both images of the pair, and manual measurements of 40 
these features confirm that their elevations vary by up to 
several htmdred meters. 

The apparent contradiction with Paillou et al.’s (2008) 
measurement has not yet been explained. One possibility is 
that the moderately-dark areas are not in fact liquid-filled seas 45 
of varying depth, but are relatively smooth, perhaps swampy, 
exposed surfaces of varying texture. This would, however, 
seem to be at odds with (a) the observation that all near- 
circular lakes appear to get progressively darker towards the 
middle until they reach the noise floor, and (b) an analysis of 50 
a channel that extends into one of the polar seas and gets 
consistently darker downstream, a result that indicates 
strongly imaging through progressively greater depths of liq- 
uid. 

The laboratory results of Paillou et al. represent only one 55 
value for absorptivity of liquid hydrocarbon that may not be 
representative of all seas and lakes, and that Titan’s seas are 
sufficiently transparent that the radar is indeed seeing the 
bottom through up to a few hundred meters of liquid. 

The seas have also been observed radiometrically. The 60 
Cassini RADAR instrument includes a passive mode able to 
acquire data with, or separately from, the active measure- 
ments. Over the past 4 years, the Cassini radiometer has 
observed the 2.2 cm thermal emission from Titan’s surface at 
resolutions ranging from 5 to 500 km, and at a variety of 65 
emission angles and polarizations. 98% of the surface has 
been mapped by the radiometer, and nearly 30% with a high 
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resolution, enabling the construction of a global map of the 
brightness temperature corrected to normal incidence (Jans- 
sen et al., 2009). 

The radar-dark seas and lakes observed at closest approach 
are well resolved on the radiometry map. They appear as high 
emissivity (radiometrically warm) regions surrounded by ter- 
rain of lower emissivity, consistent with a dielectric constant 
of about 1 .7, and an equator-to-pole temperature gradient of 
about 2 K. The brightness temperature mosaic revealed dif- 
ferences between lakes that might indicate different compo- 
sition. An accurate knowledge of the dielectric properties of 
liquid ethane/methane mixtures will allow a better interpre- 
tation of the data collected over the winter (when the evapo- 
ration can be neglected and thus the physical temperature of 
the lakes taken as constant), and may ultimately lead to a 
classification of the lakes as a function of their ethane/meth- 
ane ratio. 

Another aspect of importance is that when winter under- 
way in the northern latitudes, the volumes and mixing ratio of 
ethane and methane (the more volatile liquid) in the lakes 
might change; over the course of the Cassini mission, one 
may be observing the lakes for a total of ~1 0 terrestrial years, 
which is approximately one third of a Satum/Titan year, and 
so seasonal changes may be observed. Establishing a semi- 
empirical mixing law for the dielectric properties of a ethane/ 
methane solution will enable a better understanding of the 
impact of the seasonal changes in radar observations of the 
lakes (a major objective of the Cassini Solstice Mission). In 
addition, it will help to interpret changes in the brightness 
temperatures of the lakes by discriminating between the 
effect of a change in the composition and the effect of a 
decrease in the physical temperature. 

In addition to the atmosphere, liquids and various ices, 
Titan’s environment contains significant amounts of hydro- 
carbons, organics, tholins and possibly nitriles organic mate- 
rials resulting from photolysis of methane in the upper atmo- 
sphere, all of which may be dissolved to some extent in the 
lake liquids, and affect their dielectric properties. The present 
invention may also explore their impact on loss tangent, in 
particular, in order to assess their likely impact. 

b) Radar Modeling of Titan’s Lakes and Seas 

The methods used for determining the relationship 
between observed radar backscatter, incidence angle and lake 
floor depth, for a two-layer surface of a complex dielectric 
liquid overa solid, have already been applied by Mitchell, and 
are summarized here. In that work, the ~10 -3 loss tangent of 
Paillou et al. (2008) was used, which may be revised in view 
of the results from task (2) discussed above. 

Radar backscattered power, a 0 , is found in the uncorrected 
(for incidence angle) versions of the RADAR basic image 
data records (BIDR) products in the PDS archive. 

The present invention may fix surface and subsurface 
roughness parameters for the two-layer scattering model of 
the lake. The surface layer is characterized by its dielectric 
constant (previously taken to be£ 1 =1.75-0.002j, after lique- 
fied natural gas (LNG) measurements), surface roughness 
(height standard deviation cq, correlation length LJ, and 
albedo (a=0.2, in order to take into account scattering by 
heterogeneities or suspensions in the liquid). The bottom 
layer is characterized by its dielectric constant (E 2 =2.20- 
O.Olj for tholins after (Rodriguez et al., 2003) if non-ice 
“sediment” is suspected, and E 2 =3.10-0.000006j for water- 
ice if “bedrock” is suspected), and roughness (cr 2 , L 2 ). 

Using this method, the present invention may estimate the 
maximum observable depth on a pixel-by-pixel basis, for 
different materials in individual RADAR passes T16, T18, 
T1 9, T25, T28, T29 and T30 (all of which have been in PDS 
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archive for long enough to be eligible for OPR analysis). 
Where stereo bathymetry is available, the present invention 
may intercompare with radar modeling results for the full 
viable range of compositions. 

In doing so, the present invention may determine: (a) if 5 
these data are indeed consistent with illumination through a 
liquid body, (b) whether variations between lakes are consis- 
tent with changes in liquid composition, and (c) if possible, 
constraints by which the present invention may estimate the 
depth of observed features in non-stereo scenes. to 

During this exercise, the present invention may also look 
for possible changes of properties and depths as a function of 
latitude and elevation, drawing on 2008 Cassini Data Analy- 
sis Program (CDAP) to investigate lake elevations in the polar 
region. 15 

II. Surfaces Properties ofthe Airless Icy Satellites of Saturn 

Saturn’s airless icy satellites range from planet-like geo- 
logically active worlds to small irregular objects. Their spec- 
tra imply compositions of mostly water ice (Clark et al., 

1 986), but their respective evolutionary history and near envi- 20 
ronment have lead to different regolith composition and 
structure. Investigating the radar and radiometric properties 
of the icy satellites’ surface can help estimating their degree 
of purity and maturity as well as their possible composition, 
and therefore provide insight on the various physical and 25 
geological processes that have affected them. Yet, such effort 
requires knowledge of the microwave dielectric properties of 
the ice mixed with different contaminants at the relevant 
temperatures. 

a) Radar and Radiometric Observations of Saturnian Air- 30 
less Icy Satellites 

During the 2004 through 2007 oppositions of the Saturn 
system, Black et al. (2007) observed the mid-size Saturnian 
satellites Rhea, Dione, Tethys, and Enceladus with the 
Arecibo Observatory’s 13 cm wavelength (2.38 GHz) radar 35 
system and measured their bulk radar reflectance properties. 
Early on, Black et al. (2004) had used the same radar system 
to observe Iapetus’ optically bright, trailing hemisphere and 
the optically dark, leading hemisphere of Iapetus. 

Concurrently, since Saturn orbit insertion in 2004, the 40 
Cassini RADAR has acquired a wealth of radar data on Sat- 
urn’s main airless icy satellites, namely Mimas, Phoebe, 
Enceladus, Tethys, Dione, Rhea, Hyperion and Iapetus (Ostro 
et al., 2006; Ostro et al., 2009). With the exception of one 
close flyby of Iapetus, all these observations were distant ones 45 
and therefore designed largely for disk-integrated albedo cal- 
culations at 2.2 cm. Radiometric data obtained concurrently 
provide the thennal radiometric brightness of the icy satellites 
at the same wavelength. 

On Sep. 10, 2007, the Cassini spacecraft encountered Iape- 50 
tus at an altitude near 20,000 km and collected radar and 
radiometric measurements on the leading side of the moon, 
imaging its dark terrains and part of the bright terrains of the 
North and South polar caps. 

The datasets listed above are all publicly available, either in 55 
literature or in the PDS archive. 

b) Characterization of the Icy Satellites Surfaces and Near- 
surfaces 

At temperatures as low as 100 K, surfaces made of water 
ice, if fresh and clean, have negligible absorption length at 60 
radio wavelength (Thompson and Squyres, 1990; Maetzler, 

1 998; Lorenz, 1 998). The general understanding gained so far 
about the radar properties of the icy surfaces of the Saturn 
system, including Titan, is that they are precisely character- 
ized by relatively much lower loss compared to other solar 65 
system surfaces (Wye et al., 2007; Janssen et al., 2009). The 
dominant reflection mechanism in the microwave domain is 
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therefore expected to be subsurface volume scattering due to 
inhomogeneities within the low-loss subsurface medium. 
This mechanism has been advanced to explain the radar 
reflective properties of the Galilean satellites from Earth- 
based radar (Ostro, 1993; Black et al., 2001). 

However, the icy satellites exhibit different radar signa- 
tures and longitudinal variations that may be indicative of 
differences in the electrical properties and thus in composi- 
tion of the regolith. In particular, Ostro et al. (2006; 2009) 
invoke an increase in surface contaminants of Saturnian sat- 
ellites with distance from Saftim to explain the observed 
proportional decrease in radio albedo from Enceladus to 
Phoebe. The purity of Enceladus’ regolith may be due to its 
ongoing geologic activity and the resultant resurfacing by 
fresh water ice. 

The wavelength dependence of radar properties is a further 
argument for the likely presence of non-ice material in the icy 
regolith of the Saturnian satellites (Black et al., 2007 ; Ostro et 
al., 2006, 2009). The 13 cm wavelength radar albedos are 
systematically lower than similar results from the Cassini 
RADAR instrument at 2.2 cm, suggesting a gradient in com- 
position with depth. The most striking observation concerns 
Iapetus. This Saturnian moon is famous for the dramatic and 
unique albedo asymmetry between its leading and trailing 
hemispheres. The Cassini RADAR observations of Iapetus 
mimic the optical ones (Ostro et al., 2006), whereas the 
ground-based 1 3 cm radar measurements show hardly any 
hemispheric dichotomy (Black et al., 2004). This means that 
the dark material of Iapetus, the source of which is unknown, 
is thin. It also implies that multiple volume scattering is not as 
effective at 1 3 cm as it is at 2 cm. Black et al. (2004) and Ostro 
et al. (2009) argue that, since Iapetus’ entire surface is mature 
regolith, the wavelength dependence must involve composi- 
tion rather than structure. Either the composition is a function 
of depth everywhere (with electrical loss much greater at 
depths greater than a decimeter or two), or the regolith mate- 
rial is dispersive (with higher losses at higher wavelength). 

Both Black et al. (2004) and Ostro et al. (2009) point to 
ammonia as the most likely contaminant ofthe icy satellites’ 
regolith, because it is known to reduce the radar albedo with- 
out affecting the visual albedo (Lorenz, 1998). Ammonia is 
indeed expected to be a dominant contaminant due to the cool 
temperature in the Saturnian nebula (Lewis, 1972). It has 
been unequivocally detected on Enceladus (Emery et al., 
2005; Verbiscer et al., 2006). 

However, the present invention can think of other candi- 
dates that would increase the absorptivity of the regolith: 1) 
from cosmo-chemical considerations it can be argued that 
silicate minerals should be present on the surfaces of the icy 
satellites (Lewis, 1971), 2) organics such as nitriles, acetylene 
polymers or tholins have been detected in the spectra of a 
number of satellites (Owen et al., 2001 ; Buratti et al., 2005), 
3) the presence of metallic iron oxides, for instance, in the 
form of iron-bearing grains has been advanced (Brown et al., 
2006; Filacchione et al., 2009, submitted). The water ice can 
also be mixed with other frozen volatiles such as methane, 
carbon monoxides, nitrogen, carbon dioxide, or sulfur diox- 
ide. 

Dielectric measurements on various ice mixtures are cru- 
cial to be able to properly use the available radar and radiom- 
etry observations of the Saturn’s icy satellites to constrain the 
nature and the concentration of the contaminant in their near 
surface. These measurements will be also useful for the 
understanding of the surfaces of the Galilean satellites. 

III. Properties of Titan’s Solid Surface 

The candidates for Titan’s surface materials are similar to 
the icy satellites’ ones, though likely structurally different. 
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Plausible cryovolcanic landfonns have been identified on 
Titan (Nelson et al., 2009; Wall et al., 2009), and prompt 
measurements on ammonia-rich water ice samples (Mitri et 
al, 2008). 

However, Titan’s atmospheric and crystal volatile chemis- 5 
try produce a complex assemblage of materials, in some cases 
beyond those of other icy bodies that may be mobilized surfi- 
cially and in the upper crust. Therefore, in addition to water 
ice mixtures, candidates for Titan’s solid surface should 
include an accumulation of photochemically-produced to 
organic materials (tholins) resulting from methane in the 
upper atmosphere, and various processed materials develop- 
ing from this precipitate (as in the dunes fields). 

Resultant data may be compared to the knowledge of the 
dielectric properties of the main geological units of Titan 15 
obtained from the scatterometry analysis (Wye et al, 2007; 
Wye et al, 2008) and the radiometry polarimetric measure- 
ments (Janssen et al, 2009). 

Overview of Measurements 

Once the experimental 13.8 GHz microwave cavity is vali- 20 
dated using liquid nitrogen, whose dielectric properties are 
known (Smith et al. 1 99 1 ), the present invention may perform 
dielectric measurements on pure methane. The gas handling 
system may be consistent with safety requirements of various 
entities (e.g. JPL). 25 

From these initial measurements, the present invention 
may determine if modifications to the experimental set up are 
necessary. 

The present invention may then perform measurements on 
pure ethane and on methane/ethane mixtures covering the 30 
range of ethane from 0% to 100% in 10% intervals. 

Then, mixture measurements may be performed that 
include dissolved nitrogen and solid solutes (such as ammo- 
nia or tholins) for comparison. 

The present invention may also perform 13.8 GHz and 2.38 35 
GHz measurements on ice solids relevant to the surfaces of 
icy satellites. 

Measurements on solids may involve pure water ice (H 2 0), 
ammonia-rich water ice and water ice mixed with other fro- 
zen volatiles (methane, carbon monoxides, nitrogen, carbon 40 
dioxide, or sulfur dioxide) at the relevant temperatures. 
Samples of water ice mixed with contaminant such as sili- 
cates, organics, and iron oxides may also be studied. 

Perceived Impact to State of Knowledge 

The full interpretation of the radar and radiometry obser- 45 
vations of Saturn’s icy satellites depends critically on under- 
standing the dielectric properties of potential surface materi- 
als. Thus, the present invention may obtain and study the 
dielectric properties of relevant liquids and solids for Saturn’ s 
icy satellites. 50 

The present invention may also provide independent data 
to support (or refute) the interpretation of present-day liquid 
lakes on Titan, and may improve constrains on lake depths, 
volumes and potentially compositions, which are important 
to understanding Titan’ s carbon/organic cycle, and inevitably 55 
the evolution of its environment. 

The present invention’s measurements on solid samples 
may help constraining the surface composition of the icy 
satellites of Saturn and thus provide critical insight into the 
unique processes that have shaped them. 60 

Relevance to Element Programs and Objectives in the 
NRA 

The present invention is relevant to NASA’s strategic plan, 
to “advanced scientific knowledge of the origin and history of 
the solar system ...” and contributing to “. . . learning how the 65 
Sun’s family of planets and minor bodies originated and 
evolved.” 
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The research is relevant to the objects of “Enhancing the 
scientific return from the New Horizons, Cassini, Galileo, 
Voyager, Pioneer, and Ulysses missions by continuing the 
analysis of their respective data sets through broadened sci- 
entific participation,” “Improving our understanding of the 
evolution of the outer Solar System, including the giant plan- 
ets, their satellites, and other small bodies,” and “Creating 
data products useful to the broad scientific community.” The 
approach to be developed here promises to yield significant 
insight into Titan hydrological/hydrocarbon research. It also 
provides a clear path to future studies, specifically the analy- 
sis of the radar properties of icy worlds’ solids. 

The present invention may help to understand data from 
Cassini/Huygens studies of Titan, and other satellites that 
have only recently been released to the PDS archive. The 
present invention may help to fully appreciate the discoveries 
made by the Cassini-Huygens and Galileo missions on icy 
satellites’ surfaces, and may be of great importance to support 
the preparation of flagship missions to Europa and the Jupiter 
System, and to Enceladus and Titan. 

However, many applications are envisaged for the present 
invention. Applications for the present invention are not lim- 
ited to those applications described herein. 

Hardware Environment 

FIG. 9 is an exemplary hardware and software environment 
900 used to implement one or more embodiments of the 
invention. The hardware and software environment includes a 
computer 902 and may include peripherals. Computer 902 
may be a user/client computer, server computer, or may be a 
database computer. The computer 902 comprises a general 
purpose hardware processor 904A and/or a special purpose 
hardware processor 904B (hereinafter alternatively collec- 
tively referred to as processor 904) and a memory 906, such as 
random access memory (RAM). The computer 902 may be 
coupled to other devices, including input/output (I/O) devices 
such as a keyboard 914, a cursor control device 916 (e.g., a 
mouse, a pointing device, pen and tablet, etc.) and a printer 
928. In one or more embodiments, computer 902 may be 
coupled to a media viewing/listening device 932 (e.g., an 
MP3 player, iPod™, Nook™, portable digital video player, 
cellular device, personal digital assistant, etc.). Inoneormore 
embodiments, computer 902 may be coupled to a VNA, or 
other devices used to measure the cavity complex valued 
resonant frequencies. 

In one embodiment, the computer 902 operates by the 
general purpose processor 904A performing instructions 
defined by the computer program 910 under control of an 
operating system 908. The computer program 910 and/or the 
operating system 908 may be stored in the memory 906 and 
may interface with the user and/or other devices to accept 
input and commands and, based on such input and commands 
and the instructions defined by the computer program 910 and 
operating system (OS) 908 to provide output and results. 

Output/results may be presented on the display 922 or 
provided to another device for presentation or further pro- 
cessing or action. In one embodiment, the display 922 com- 
prises a liquid crystal display (LCD) having a plurality of 
separately addressable liquid crystals. Each liquid crystal of 
the display 922 changes to an opaque or translucent state to 
form a part of the image on the display in response to the data 
or information generated by the processor 904 from the appli- 
cation of the instructions of the computer program 910 and/or 
operating system 908 to the input and co mm ands. The image 
may be provided through a graphical user interface (GUI) 
module 918A. Although the GUI module 918A is depicted as 
a separate module, the instructions performing the GUI func- 
tions can be resident or distributed in the operating system 
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908, the computer program 910, or implemented with special 
purpose memory and processors. 

Some or all of the operations performed by the computer 
902 according to the computer program 910 instructions may 
be implemented in a special purpose processor 904B. In this 5 
embodiment, the some or all of the computer program 910 
instructions may be implemented via firmware instructions 
stored in a read only memory (ROM), a programmable read 
only memory (PROM) or flash memory within the special 
purpose processor 904B or in memory 906. The special pur- to 
pose processor 904B may also be hardwired through circuit 
design to perform some or all of the operations to implement 
the present invention. Further, the special purpose processor 
904B may be a hybrid processor, which includes dedicated 
circuitry for performing a subset of functions, and other cir- 15 
cuits for performing more general functions such as respond- 
ing to computer program instructions. In one embodiment, 
the special purpose processor is an application specific inte- 
grated circuit (ASIC). 

The computer 902 may also implement a compiler 912 20 
which allows an application program 910 written in a pro- 
gramming language such as COBOL, Pascal, C++, FOR- 
TRAN, or other language to be translated into processor 904 
readable code. After completion, the application or computer 
program 910 accesses and manipulates data accepted from 25 
I/O devices and stored in the memory 906 of the computer 
902 using the relationships and logic that was generated using 
the compiler 912. 

The computer 902 also optionally comprises an external 
communication device such as a modem, satellite link, Eth- 30 
emet card, or other device for accepting input from and pro- 
viding output to other computers 902. 

In one embodiment, instructions implementing the operat- 
ing system 908, the computer program 910, and the compiler 
912 are tangibly embodied in a computer-readable medium, 35 
e.g., data storage device 920, which could include one or 
more fixed or removable data storage devices, such as a zip 
drive, floppy disc drive 924, hard drive, CD-ROM drive, tape 
drive, etc. Further, the operating system 908 and the computer 
program 91 0 are comprised of computer program instructions 40 
which, when accessed, read and executed by the computer 
902, causes the computer 902 to perform the steps necessary 
to implement and/or use the present invention or to load the 
program of instructions into a memory, thus creating a special 
purpose data structure causing the computer to operate as a 45 
specially programmed computer executing the method steps 
described herein. Computer program 910 and/or operating 
instructions may also be tangibly embodied in memory 906 
and/or data communications devices 930, thereby making a 
computer program product or article of manufacture accord- 50 
ing to the invention. As such, the terms “article of manufac- 
ture,” “program storage device” and “computer program 
product” as used herein are intended to encompass a com- 
puter program accessible from any computer readable device 
or media. 55 

Of course, those skilled in the art will recognize that any 
combination of the above components, or any number of 
different components, peripherals, and other devices, may be 
used with the computer 902. 

Although the term “user computer” or “client computer” is 60 
referred to herein, it is understood that a user computer 902 
may include portable devices such as cell phones, notebook 
computers, pocket computers, or any other device with suit- 
able processing, communication, and input/output capability. 

FIG. 10 schematically illustrates a typical distributed com- 65 
puter system 1000 using a network 1002 to connect client 
computers 902 to server computers 1006. A typical combina- 
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tion of resources may include a network 1002 comprising the 
Internet, LANs (local area networks), WANs (wide area net- 
works), SNA (systems network architecture) networks, or the 
like, clients 902 that are personal computers or workstations, 
and servers 1006 that are personal computers, workstations, 
minicomputers, or mainframes (as set forth in FIG. 9). 

A network 1002 such as the Internet connects clients 902 to 
server computers 1006. Network 1002 may utilize ethemet, 
coaxial cable, wireless communications, radio frequency 
(RF), etc. to connect and provide the communication between 
clients 902 and servers 1006. Clients 902 may execute a client 
application or web browser and communicate with server 
computers 1006 executing web servers 1010. Such a web 
browser is typically a program such as MICROSOFT 
INTERNET EXPLORER™, MOZILLA FIREFOX™, 
OPERA™, APPLE SAFARI™, etc. Further, the software 
executing on clients 902 may be downloaded from server 
computer 1006 to client computers 1002 and installed as a 
plug in or ACTIVEX™ control of a web browser. Accord- 
ingly, clients 902 may utilize ACTIVEX™ components/com- 
ponent object model (COM) or distributed COM (DCOM) 
components to provide a user interface on a display of client 
902. The web server 910 is typically a program such as 
MICROSOFT'S INTERNET INFORMATION SERVER™. 

Web server 1010 may host an Active Server Page (ASP) or 
Internet Server Application Programming Interface (ISAPI) 
application 1012, which may be executing scripts. The scripts 
invoke objects that execute business logic (referred to as 
business objects). The business objects then manipulate data 
in database 1016 through a database management system 
(DBMS) 1014. Alternatively, database 1016 may be part of or 
connected directly to client 902 instead of communicating/ 
obtaining the information front database 1016 across network 
1002. When a developer encapsulates the business function- 
ality into objects, the system may be referred to as a compo- 
nent object model (COM) system. Accordingly, the scripts 
executing on web server 1010 (and/or application 1012) 
invoke COM objects that implement the business logic. Fur- 
ther, server 1006 may utilize MICROSOFT’S™ Transaction 
Server (MTS) to access required data stored in database 1016 
via an interface such as ADO (Active Data Objects), OLE DB 
(Object Linking and Embedding DataBase), or ODBC (Open 
DataBase Connectivity). 

Generally, these components 1 006-1 01 6 all comprise logic 
and/or data that is embodied in/or retrievable from device, 
medium, signal, or carrier, e.g., a data storage device, a data 
communications device, a remote computer or device 
coupled to the computer via a network or via another data 
communications device, etc. Moreover, this logic and/or data, 
when read, executed, and/or interpreted, results in the steps 
necessary to implement and/or use the present invention 
being performed. 

Although the term “user computer”, “client computer”, 
and/or “server computer” is referred to herein, it is under- 
stood that such computers 902 and 1006 may include portable 
devices such as cell phones, notebook computers, pocket 
computers, or any other device with suitable processing, com- 
munication, and input/output capability. 

Of course, those skilled in the art will recognize that any 
combination of the above components, or any number of 
different components, peripherals, and other devices, may be 
used with computers 902 and 1006. 

Software Embodiments 

Embodiments of the invention are implemented as a soft- 
ware application on a client 902 or server computer 1006. 
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Conclusion 

35 This concludes the description of the preferred embodi- 
ment of the present invention. The foregoing description of 
one or more embodiments of the invention has been presented 
for the purposes of illustration and description. It is not 
intended to be exhaustive or to limit the invention to the 
40 precise form disclosed. Many modifications and variations 
are possible in light of the above teaching. It is intended that 
the scope of the invention be limited not by this detailed 
description, but rather by the claims appended hereto. 

What is claimed is: 

45 1 . A method of calculating, in one or more computers, a 

complex-valued dielectric permittivity (£) or a complex val- 
ued magnetic permeability (p) of at least one material, com- 
prising: 

(a) receiving, in one or more of the computers, a measured 
50 complex-valued resonance frequency Fmaterial of a 

cavity containing the material, wherein Fmaterial is a 
measurement obtained by exciting the cavity with one or 
more transverse magnetic TM 0 „ 0 cavity resonant modes ; 

(b) in one or more of the computers: 

55 (1) receiving a solution of Maxwell’s equations for a 

complex-valued resonance frequency function F of 
the cavity containing the material, wherein the F is a 
function of £ or p or £ and p of the material, 

(2) numerically or analytically solving an equation 
60 F=Fmaterial exactly for the £ or the p, wherein the £ 

or the p is the only unknown; and 

(c) outputting, from one of the computers, the £ or p. 

2 . The method of claim 1 , wherein the cavity and the 
material are annular or cylindrical regions. 

65 3. Hie method of claim 2, wherein the material is non- 

magnetic, has any cylindrical size, and is positioned along the 
cavity’s axis. 
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4. The method of claim 2, wherein to obtain the p of the 
material that is magnetic, the cavity resonant modes comprise 
at least one Transverse Magnetic (TM 0 „ 0 ) mode, and a mag- 
netic field and an electric field associated with the TM 0 „ 0 
mode, where n is an integer, the method further comprising: 

(i) performing the steps (a)-(b) for the material that is a first 
cylindrical material positioned along the cavity’s axis, 
the first cylindrical material having a first radius that is 
sufficiently small so that the magnetic field’s strength at 
the first radius’s end is negligible, thereby allowing the 
Maxwell’s equations and the F to be solved for only one 
unknown, the E=Ematerial; and 

(ii) performing the steps (a)-(b) for the material that is a 
second cylinder comprising a same material as the first 
cylinder, and having a larger radius than the first radius, 
such that the second cylinder positioned along the axis 
experiences both the magnetic field and the electric field, 
thereby allowing the Maxwell’s equations and the F to 
be solved for only one unknown, the p, using the 
E=Ematerial obtained in step (i). 

5. The method of claim 2, wherein to obtain the p of the 
material that is magnetic, the cavity resonant modes com- 
prise: 

a first frequency, a first magnetic field, and a first electric 
field associated with a first TM 0 „ 0 mode, where n is an 
integer, 

a second frequency, a second magnetic field and a second 
electric field associated with a second TM 0fl0 mode, 
the first frequency and the second frequency such that the E 
and the p at the first frequency are negligibly different 
from the E and the p at the second frequency; and 
the method further comprising: 

(i) exciting the cavity with the first TM 0 „ 0 mode where the 
material is an annulus positioned at a radius of the cavity 
associated with a maximum of the first electric field and 
a zero of the first magnetic field; 

(ii) performing the steps (a)-(b) using the first TM 0w0 mode, 
thereby allowing the Maxwell’s equations and the F to 
be solved for only one unknown, the Ematerial; and 

(iii) performing steps (a)-(b) using the second TM 0w0 mode 
so that the material experiences both the second electric 
field and the second magnetic field, thereby allowing the 
Maxwell’s equations and the F to be solved for only one 
unknown, the p, using the E=Ematerial obtained in step 
(ii). 

6. The method of claim 2, further comprising: 
receiving, in one or more of the computers, a measured 

complex-valued resonance frequency Fcavity of the 
cavity without the material; 

in one or more of the computers, modifying the F to include 
a cavity wall loss contribution associated with a change 
in the transverse magnetic cavity modes’ fields at the 
cavity’s walls and electrical conductivity ct„, of the cavi- 
ty’s walls, calculated from Fcavity and Fmaterial; and 
in the step (b)(2), numerically or analytically solving the 
equation F=Fmaterial, wherein the F is the modified 
function F. 

7. The method of claim 6, further comprising: 
calculating, in one or more of the computers: 

the effective electrical conductivity a w of the cavity’s 
walls from 

a w | Q 0 2 /(jtp 0 f 0 ) |(2=D ( /L) 2 /I ), 2 , where Q 0 =Re (Fcav- 
ity )/2 Im(Fcavity), p,, is the permeability of vacuum, 
f 0 =Re(Fcavity), D c is the cavity diameter, and L is the 
cavity length, wherein the Fcavity is the measured 
complex-valued resonance frequency of the cavity 
without the material; 
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R e =[(3tp 0 Re[Fmaterial])A. w ] 1/2 ; 

P M , a/? =(R J /2)/ s IFl / l 2 ds, where indicates an integral over 
the surface of the cavity walls, ds is the surface ele- 
ment, and FI r is the magnetic field tangential compo- 

5 nent at the surface of the cavity walls for the cavity 

containing the material; 

Q wa „=2jt Re(Fmaterial) U/P wa/; where U is the total 
electromagnetic energy stored in the cavity contain- 
ing the material; and 

to the cavity wall loss contribution f w; =Re (Fmaterial)/2 
Q and 

modifying the F into F=F+jf vra// to account for the power 
P waH absorbed by the cavity’s walls associated with the 
electrical conductivity; and 

15 solving the equation F=F +jf vra/7 =Fmaterial for the E or the 
p in the step b(2). 

8 . The method of claim 6, wherein the material in the cavity 
is in a material holder, and the material holder has a complex- 
valued dielectric permittivity (Eholder), the method further 

20 comprising: 

obtaining a measured complex-valued resonance fre- 
quency Fholder of the cavity comprising the material 
holder without without any material in the material 
holder; 

25 solving the Maxwell’s equations for a complex-valued 
resonance frequency g of the cavity including the mate- 
rial holder without any material in the material holder, 
wherein the g is a function of the Eholder; 

modifying the g to account for the cavity’s wall losses 

30 associated with the ow, thereby obtaining a resultant 
complex -valued resonant frequency G that is a function 
of the Eholder; 

obtaining the Eholder by solving G=Fholder; and 

performing the steps (a)-(b) for the material inserted in the 

35 material holder in the cavity, wherein the solving of 
F=Fmaterial uses the Eholder obtained by the solving of 
G=Fholder. 

9. The method of claim 8, further comprising measuring 
the Fholder, Fcavity, and Fmaterial wherein the material is a 

40 liquid, powder, or granular material, 

10. The method of claim 6, further comprising measuring 
the Fcavity and the Fmaterial. 

11 . The method of claim 1 , wherein the solving further 
comprises: 

45 (i) obtaining or selecting one or more values for the E or the 

p; 

(ii) calculating the F for one or more of the values; 

(iii) comparing the F with Fmaterial; 

repeating steps (i)-(iii) with another of the values until 

50 the value that satisfies F=Fmaterial to any desired precision 
is selected. 

12 . The method of claim 1 , wherein the cavity comprises 
one or more additional materials having known material 
dielectric permittivity and known material magnetic penne- 

55 ability, such that the E or p is an only unknown during the 
solving step (b). 

13 . The method of claim 1 , wherein the solving of 
F=Fmaterial further comprises finding roots of F-Fmate- 
rial=0, by: 

60 (i) specifying a range of values for z within a first rectangle 

in a complex plane, wherein the z are complex numbers 
in a complex plane; 

(ii) introducing a mesh into the first rectangle, thereby 
forming the mesh comprising mesh points at corners of 

65 second rectangles, wherein the second rectangles are 
smaller than the first rectangles; 

(iii) evaluating Re F(z) at each of the mesh points; 
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(iv) determining one or more first sides of the second 
rectangles where the Re (F(z)) changes sign; 

(v) evaluating Im F(z) at each of the mesh points; 

(vi) determining one or more second sides of the second 
rectangles where the Im (F(E)) changes sign; 

(vii) selecting the second rectangles having both the first 
sides and second sides, where both the Re F(z) and the 
Im F(z) change sign, to obtain selected second rect- 
angles; 

(viii) interpolating the Re F(z) in the selected second rect- 
angles to define a Re F(z) null line; 

(ix) interpolating the Im F(z) in the selected second rect- 
angles to define an Im F(z) null line; 

(x) finding a first intersection point between the Re F(z) 
null line and the Im F(z) null line; 

(xi) setting the second rectangle as the first rectangle and 
repeating steps (i)-(x) using the second rectangle as the 
first rectangle, thereby finding a second intersection 
point; 

(xii) obtaining a difference between the first intersection 
point and the second intersection point, and comparing 
the difference with a precision parameter that is equal to 
a desired precision; 

(xiii) repeating the method until the difference is at least as 
small as the precision parameter, thereby obtaining a 
final intersection point; and 

(xiv) setting the final intersection point equal to the E or p. 

14. The method of claim 1, further comprising: 

using the E or p to interpret remote sensing data obtained 
from earth or extra-terrestrial material. 

15. A non-transitory computer readable storage medium 
encoded with computer program instructions which when 
accessed by a computer cause the computer to load the pro- 
gram instructions to a memory therein creating a special 
purpose data structure causing the computer to operate as a 
specially programmed computer, executing a method of cal- 
culating a complex-valued dielectric permittivity (E) or a 
complex valued magnetic permeability (p) of at least one 
material, comprising: 

(a) receiving, in the specially programmed computer, a 
measured complex-valued resonance frequency Fcavity 
of a cavity without the material; 

(b) receiving, in the specially programmed computer, a 
measured complex -valued resonance frequency Fmate- 
rial of the cavity containing the material, wherein Fcav- 
ity and Fmaterial are measurements obtained by exciting 
the cavity with one or more transverse magnetic TM 0 „ 0 
cavity resonant modes; and 

(c) numerically or analytically solving, in the specially 
programmed computer, Maxwell’s equations exactly for 
E or p wherein the solving includes: 

(i) numerically or analytically calculating a complex 
valued resonance frequency F for each of one or more 
values of E or p, wherein the F includes a cavity wall 
loss contribution associated with a change in the tran- 
verse magnetic cavity modes’ fields at cavity’s walls 
and calculated from Fcavity and Fmaterial; and 

(ii) selected the E or p wherein F=Fmaterial, therby 
solving for the E or p; and 

(d) outputting, from the specially programmed computer, 
the selected E or p. 

16. The computer readable storage medium of claim 15, 
wherein the cavity and the material are annular or cylindrical 
regions. 
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17. The computer readable storage medium of claim 15, 
wherein, the calculating in step c(i) further comprises: 

calculating, in the specially programmed computer: 
an effective electrical conductivity c M , of the cavity’s 
5 walls from 

a v ,=[Qo 2 /Uhofo)](2=D^L) 2 /D <: 2 , where Q 0 =Re (Fcav- 
ity )/2 Im(Fcavity), p,, is the permeability of vacuum, 
f 0 =Re(Fcavity), is the cavity diameter, and L is the 
cavity length; 

to R e =[(jtp 0 Re| Fmaterial])/^,] 1/2 ; 

Pv,.a// = (R J /2)/ s IH / l 2 ds, where indicates an integral over 
the surface of the cavity walls, ds is the surface ele- 
ment, and H, is the magnetic field tangential compo- 
nent at the surface of the cavity walls for the cavity 
15 containing the material; 

Q w/ =2jTRe(Fmaterial) U/P M , a „ where U is the total elec- 
tromagnetic energy stored in the cavity containing the 
material; and 

the cavity wall loss contribution f waH =Re (Fmaterial)/2 
20 Q wall to account for the P M , a// power absorbed by the 

cavity’s walls; and 

the solving in the step c(ii) comprises solving F=F+ 

jU<i// = F mater ial- 

18. The computer readable storage medium 17, wherein 
25 the material in the cavity is in a material holder, and the 

material holder has a complex-valued dielectric permittivity 
(Eholder), the method further comprising: 

receiving, in the specially programmed computer, a mea- 
sured complex-valued resonance frequency Fholder of 
30 the cavity using a TM 0 „ 0 mode, the cavity comprising 
the material holder without any material in the material 
holder; 

solving, in the specially programmed computer, the Max- 
well’s equations for a complex-valued resonance fre- 
35 quency g of the cavity including the material holder 
without any material in the material holder, wherein the 
g is a function of the Eholder; 

modifying, in the specially programmed computer, the g to 
account for the electrical conductivity, thereby obtaining 
40 a resultant complex -valued resonant frequency G that is 
a fimction of the Eholder; 

obtaining, in the specially programmed computer, the 
Eholder by solving G=Fholder; and 

performing the steps (a)-(c) for the material inserted in the 
45 material holder in the cavity, wherein the solving of 
F=Fmaterial uses the Eholder obtained by the solving of 
G=Fholder. 

19. The computer readable storage medium of claim 15, 
fiirther comprising: 

50 receiving, in the specially programmed computer, a range 
of values for the E or p; and 

selecting, in the specially programmed computer, the value 
in the range that satisfies F=Fmaterial to any degree of 
precision. 

55 20. The computer readable storage medium of claim 15, 

wherein the material is non-magnetic, has any cylindrical 
size, and is positioned along an axis of the cavity. 

21. The computer readable storage medium of claim 15, 
wherein to obtain the p of the material that is magnetic, the 
60 cavity resonant modes comprise at least one Transverse Mag- 
netic (TM 0n0 ) mode, and a magnetic field and an electric field 
associated with the TM 0;;0 mode, where n is an integer, the 
method fiirther comprising: 

( 1 ) performing the steps (a)-(c ) for the material that is a first 
65 cylindrical material positioned along an axis of the cav- 
ity, the first cylindrical material having a first radius that 
is sufficiently small so that the magnetic field’s strength 
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at the first radius’s end is negligible, thereby allowing 
the Maxwell’s equations and the F to be solved for only 
one unknown, the E; and 

(2) performing the steps (a)-(c) for the material that is a 
second cylinder comprising a same material as the first 5 
cylinder, and having a larger radius than the first radius, 
such that the second cylinder positioned along the axis 
experiences both the magnetic field and the electric field, 
thereby allowing the Maxwell’s equations and the F to 
be solved for only one unknown, the p, using the E to 
obtained in step (1). 

22. The computer readable storage medium of claim 15, 
wherein to obtain the p of the material that is magnetic, the 
cavity resonant modes comprise: 

a first frequency, a first magnetic field, and a first electric 15 
field associated with a first TM 0jl0 mode, where n is an 
integer, 

a second frequency, a second magnetic field and a second 
electric field associated with a second TM 0 „ 0 mode, 

the first frequency and the second frequency such that the E 20 
and the p at the first frequency are negligibly different 
from the E and the p at the second frequency; and 

the method further comprising: 

(1 ) exciting the cavity with the first TM 0w0 mode where the 
material is an annulus positioned at a radius of the cavity 25 
associated with a maximum of the first electric field and 
a zero of the first magnetic field; 

(2) performing the steps (a)-(c) using the first TM 0 , ;0 mode, 
thereby allowing the Maxwell’s equations and the F to 
be solved for only one unknown, the E; and 

(3) performing steps (a)-(c) using the second TM 0 „ 0 mode 
so that the material experiences both the second electric 
field and the second magnetic field, thereby allowing the 
Maxwell’s equations and the F to be solved for only one 
unknown, the p, using the E obtained in step (2). 

23. The computer readable storage medium of claim 15, 
wherein the cavity comprises one or more additional materi- 
als having known material dielectric permittivity and known 
material magnetic permeability, such that the E or p is an only 
unknown during the solving step (c). 

24. The computer readable storage medium 15, wherein 
the solving of F=Fmaterial further comprises finding roots of 
F-Fmaterial=0, by: 

(i) specifying, in the specially programmed computer, a 
range of values for z within a first rectangle in a complex 45 
plane, wherein the z are complex numbers in a complex 
plane; 

(ii) introducing, in the specially programmed computer, a 
mesh into the first rectangle, thereby forming the mesh 
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comprising mesh points at corners of second rectangles, 
wherein the second rectangles are smaller than the first 
rectangles; 

(iii) evaluating, in the specially programmed computer, Re 
F(z) at each of the mesh points; 

(iv) determining, in the specially programmed computer, 
one or more first sides of the second rectangles where the 
Re F(z) changes sign; 

(v) evaluating, in the specially programmed computer, Im 
F(z) at each of the mesh points; 

(vi) determining, in the specially programmed computer, 
one or more second sides of the second rectangles where 
the Im F(E) changes sign; 

(vii) selecting, in the specially programmed computer, the 
second rectangles having both the first sides and second 
sides, where both the Re F(z) and the Im F(z) change 
sign, to obtain selected second rectangles; 

(viii) interpolating, in the specially programmed computer, 
the Re F(z) in the selected second rectangles to define a 
Re F(z) null line; 

(ix) interpolating, in the specially programmed computer, 
the Im F(z) in the selected second rectangles to define an 
Im F(z) null line; 

(x) finding, in the specially programmed computer, a first 
intersection point between the Re F(z) null line and the 
Im F(z) null line; 

(xi) setting, in the specially programmed computer, the 
second rectangle as the first rectangle and repeating 
steps (i)-(x) using the second rectangle as the first rect- 
angle, thereby finding a second intersection point; 

(xii) obtaining, in the specially programmed computer, a 
difference between the first intersection point and the 
second intersection point, and comparing the difference 
with a precision parameter that is equal to a desired 
precision; 

(xiii) repeating, in the specially programmed computer, the 
method until the difference is at least as small as the 
precision parameter, thereby obtaining a final intersec- 
tion point; and 

(xiv) setting, in the specially programmed computer, the 
final intersection point equal to the E or |x. 

25. The computer readable storage medium of claim 15, 
further comprising: 

using, in the specially programmed computer, the E or (X to 
interpret remote sensing data obtained from earth or 
extra-terrestrial material. 



